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CONFIGURATION MANAGEMENT AND AUTOMATIC CONTROL OF AN AUGMENTOR 


WING AIRCRAFT WITH VECTORED THRUST 
Luigi S. Cicolani, B. Sridhar,* and George Meyer 
Ames Research Center 


SUMMARY 


An advanced structure for automatic flight-control logic for powered-lift 
aircraft operating in terminal areas is under investigation at Ames Research 
Center. This structure is based on acceleration control; acceleration commands 
are constructed as the sum of acceleration on the reference trajectory and a 
corrective feedback acceleration to regulate path tracking errors. The central 
element of the structure, termed a Trimmap, uses a model of the aircraft 
aerodynamic and engine forces to calculate the control settings required to 
generate the acceleration commands. The feedback control design for nonlinear 
plants is simplified by this structure, since the transfer function of the 
controlled element (the Trimmap-aircraf t combination) is approximately unity 
at all configurations and flight conditions. Further, consideration of plant 
nonlinearity is isolated in the Trimmap and the nonlinearity is algebraic 
rather than dynamic. 

This report describes the design criteria for the Trimmap and derives a 
Trimmap for Ames' experimental augmentor wing jet STOL research aircraft. The 
principal problems are associated with control redundancy (there are two more 
controls than necessary to generate any given acceleration command) and model 
nonlinearity. Control redundancy is resolved using a stored configuration 
schedule which selects two of the controls (flap and engine exhaust nozzle) 
as a function of the reference flight condition, while the remaining two 
controls (elevator and throttle) are computed on the basis of the total 
acceleration command and the aircraft force model. The configuration 
schedule is designed off-line to optimize configuration-dependent parameters 
of interest subject to operational constraints on controls usage; the present 
concept maximizes lift margin and regulator control margin. The algebraic 
nonlinearity of the force model is treated using piecewise linear 
descriptions of the model over its finite domain; thus, model accuracy is 
controlled by the number of pieces and the model is inverted by solving 
linear equations for each of a finite number of pieces. 

The automatic control system, including the Trimmap described in this 
paper, was subjected to simulation tests using a rigorous STOL approach 
trajectory. These tests demonstrate the system response to low-frequency 
inputs, including maneuver commands, Trimmap model errors, and steady winds. 

The proposed Trimmap maintains an optimal configuration and coordinates all 
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controls during any maneuvers within the terminal area operational domain 
and in response to steady winds; and it compensates automatically for model 
errors. 


INTRODUCTION 


An effort is under way at Ames Research Center to derive methodology for 
automatic flight-control logic for powered-lift V/STOL aircraft flying in the 
terminal approach area and to apply the methodology to specific aircraft to 
demonstrate its feasibility. 

Among the powered-lift concepts currently being developed are the 
augmentor wing, tilt rotor, lift fan, and externally blown flap. When com- 
pared to conventional aircraft, these aircraft pose several special flight- 
control logic problems. First, the use of powered lift at landing entails 
large changes in the configuration of the aircraft from its conventional 
cruise configuration to the powered-lift landing configuration. As a result, 
aerodynamic forces are configuration- and power-dependent and are sufficiently 
nonlinear that the conventional method of approximating a nonlinear plant by 
a series of linear perturbation models at different nominal flight conditions 
and configurations may require prohibitively large numbers of such models in 
order to provide stable, near-optimum transient response over the full flight 
envelope. Even at a single flight condition, the variation of forces with 
controls may be sufficiently nonlinear that nonlinear transient response is 
obtained from linear control laws. In addition, V/STOL aircraft achieve 
powered lift by a wide variety of design concepts with radically different 
control characteristics and aerodynamic force models; consequently, the 
methodology should accommodate these variations systematically. 

A second complicating feature is control redundancy; that is, some STOL 
aircraft have a novel configuration control such as the engine exhaust nozzles 
of Ames’ augmentor wing jet STOL research aircraft (AWJSRA) or the mast angle 
of the tilt rotor. Consequently, there are many control settings which trim 
the aircraft at a given flight condition and, in fact, for the AWJSRA the set 
of such trim settings has two degrees of freedom. The automatic control 
system must resolve this redundancy and do so automatically and optimally for 
all paths within the approach flight regime. 

Third, the response of aircraft controls will vary considerably with 
configuration so that the use of these controls to regulate path errors must 
also vary with configuration during an approach. For redundantly controlled 
aircraft the novel control may also be suitable for regulator usage and 
thereby provide additional control margin. 

Finally, it is anticipated that STOL aircraft will operate in terminal 
areas in conjunction with advanced air traffic control (ATC) and guidance 
systems which make full use of STOL capability for steep descents and low- 
speed maneuverability and will command four -dimensional (time of arrival and 
path) approaches in limited airspace. The output set of these systems is 
assumed to consist of any trajectory within the flight envelope limits of the 
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aircraft, and within passenger comfort and operational constraints. The 
automatic flight-control system should be able to couple the powered-lift STOL 
aircraft to such a system and the design methodology must accommodate its full 
output set. In view of the complexities involved in flying these aircraft, 
either manually or automatically, it can be assumed that an airborne digital 
computer and associated control actuators will be a key element in V/STOL 
flight control and the methodology can assume such a digital system. 


Automatic Control System Design Structure 

A logical structure to accommodate these problems is proposed in refer- 
ence 1. A block diagram outlining the main elements in that structure for 
its application to the AWJSRA is shown in figure 1. The system is to be 
capable of executing four-dimensional trajectories as required by an advanced 
ATC system controlling terminal area STOL operations, and these may be any 
trajectories within a domain bounded by aircraft performance limits, passenger 
comfort constraints, and other operational considerations. These are assumed 
to be coarsely defined commands; for example, sequences of straight line and 
circular-arc legs. The command generator element of the logic is an interface 
whose function is to modify the ATC command as necessary to satisfy 
operating constraints on maneuvering and on the limits of aircraft maneuver 
capabilities, and which minimizes the maneuvering activity required to follow 
the ATC command. Theoretical development of the command generator and 
simulation test results are presented in reference 2. 

The output of the command generator is a time-varying nominal trajectory: 

|(RCj.(t), VCj.(t) ,aCj.(t)) , te [tQ,tp]} 

which is assumed given in inertial or runway coordinates. ^ The remaining 
elements of the control logic — the path regulator. Force Trimmap, and 
attitude control — track these commands. A special feature of the control 
structure shown in figure 1 is that the path is tracked on the basis of 
acceleration commands. The nominal path acceleration, ac^ , and a corrective 
feedback acceleration, Aa^, calculated from path errors, are summed to form 
the total acceleration command, act^. The regulator control law, 

Aa^(RCr - R^, VC^ - , operates on errors from the reference path; it is 

linear for small errors with limiting in accordance with the desired accelera- 
tion authority for the regulator and acceleration rate capabilities of the 
aircraft. This control law specifies the desired transient dynamics of the 
aircraft and control-system combination and these dynamics are achieved in 
accordance with the accuracy with which the remaining control system tracks 
the acceleration commands. 


^The reference frame in which a vector is expressed is indicated by 
subscripts; the subscripts r, b, s, and p refer, respectively, to runway, 
body, stability, and path axes. When a reference frame is specified, the 
superscript bar is dropped from the vector notation. 
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The control settings necessary to generate the commanded acceleration, 
ac^, are calculated in the Force Trimmap element of figure 1 as solutions 
of the aircraft force-balance or "trim** equations and additional kinematic 
relations. The Trimmap outputs are commands directly to control servos 
{6f, 6|-, V, and to an attitude control subsystem {(|), 6, i(j, These 

latter are Euler attitude angles and body axis components of angular velocity. 
Denoting a model of the aircraft engine and aerodynamic forces as: 

a = f(v,u,p) (1) 

where f neglects plant dynamics and is based on available data from theory, 
wind-tunnel tests or flight tests, then the Trimmap is a partial inverse of 
this model and has the form: 


= g(ac,v,p) (2) 

The Trimmap provides dynamic trimming of the aircraft because ac^ is time- 
varying during an approach, but it can be assumed that the bandwidth of ac^ 
is low compared to the bandwidth of the force dynamics and that this is 
imposed in the system by appropriate desirgn of the command generator and regu- 
lator elements of the structure. 

The Trimmap outputs include the attitude states required by the trajectory 
and these are executed by the attitude control subsystem of figure 1 using 
elevator, wheel, and rudder commands. Although this subsystem is not outlined 
in detail in figure 1, its proposed structure in reference 1 parallels that 
of the outer loop, with an attitude command generator to modify coarse 
attitude commands from the Trimmap, and an attitude regulator and moment 
Trimmap to track these commands with the desired attitude transient response 
by commanding and controlling angular acceleration. 

Finally, a state estimation system completes the control loop by provid- 
ing estimates (R^, V^, a^, G, 03^) of aircraft dynamic response 

to the servo commands and wind disturbances. 

The proposed automatic control structure is suited to powered-lift 
aircraft since the element controlled by the outer loop is the combined Force 
Trimmap (eq.. (2)) and aircraft acceleration response (eq. (1), approximately) 
which has a transfer function of approximately unity everywhere in the flight 
envelope at sufficiently low input frequencies. In that case, the regulator 
control law gains are fixed and the aircraft response is linear and invariant 
with flight condition and configuration. The nonlinearity in the control 
structure is that of the Trimmap (eq. (2)), but this is now an algebraic 
nonlinearity rather than a dynamic one. Analogous properties occur for the 
attitude control subsystem of reference 1 since the element controlled by the 
inner loop is the combined moment Trimmap and aircraft angular acceleration 
response, and this also has a transfer function of approximately unity. 

A Trimmap for the AWJSRA is derived in this report and some simulation 
results are presented. Although the Trimmap is specialized to each aircraft, 
especially in the case of V/STOL aircraft, an effort has been made to formulate 
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concepts, problem statements, and solution methods for general applicability 
to powered-lift, redundantly-controlled aircraft. 

A flow chart of the Trimmap computations for the AWJSRA is shown in fig- 
ure 2. The runway axis-referenced acceleration command, act^, is combined 
with gravity to obtain the commanded applied specific force and this is 
mapped into stability axis components (A^j, A]^)^, along and normal to the 
flight path, and the roll angle, (1)^, which locates the normal force relative 
to the vertical plane. This mapping requires the velocity vector direction 
angles (ya> ^VA^ sideslip command, 3c > if nonzero. The specific force 

commands (A^^^, Aj^^) are inputs to the Basic Trimmap element of the diagram, 
which solves the stability axis force balance equations of the aircraft. 

These are the familiar lift and drag trim equations and, together with addi- 
tional relationships to resolve the control redundancy, provide solutions for 
the required angle of attack, Uq, and the control commands (6fc> ^c> ^tc^ • 

The remaining Trimmap outputs are the trim elevator and attitude commands. 

The trim elevator is that setting which nulls the pitch angular acceleration 
and is a function of the output variables of the basic Trimmap. In the 
attitude control subsystem, the pitch axis is controlled using variations 
of the elevator about this setting. The required Euler attitude angles 
($Cj computed from angles associated with the trajectory velocity 

and acceleration and with the trim solution (ya> ^VA» ^ c ) while the 

angular velocity command is computed approximately from the rotations of the 
commanded velocity vector (y, • (More generally, the derivatives of 

(Ya> ^VA’ ^v» ^c» ^c) t)e used to construct the required angular velocity 

and angular acceleration vectors.) 

In general, the Force Trimmap for any STOL or CTOL aircraft can be 
represented as in figure 2 by changing the control variable notation as neces- 
sary. The mapping to stability axes commands and the computation of attitude 
commands use general kinematic relations that are aircraft-independent. 
Aircraft force model details appear only in the basic Trimmap; in the case of 
CTOL aircraft this model is relatively simple and the corresponding basic 
Trimmap is readily derived. For powered-lift, redundantly-controlled 
aircraft, model complexity is the source of the principal problems associated 
with the Force Trimmap; in the AWJSRA example of this report these stem 
principally from control redundancy. The proposed structure for dealing with 
this redundancy is shown in the basic Trimmap of figure 2 and consists of two 
elements: an automatic configuration schedule (ACS) and a trim solution 

algorithm. The configuration schedule is a stored schedule for two configura- 
tion controls which is calculated off-line to optimize the configuration 
point-by-point over the operational flight envelope subject to operational 
constraints on controls usage. In flight, the configuration is commanded 
open-loop by the ACS as a function of the commanded flight condition and thus 
provides optimized configuration commands along any terminal area trajectory 
commanded by the command generator. For the AWJSRA, the two controls that 
are scheduled are flap and nozzle; the flap is selected to maximize lift 
margin and the nozzle is selected to maximize the control margin of the 
remaining two controls, elevator and throttle, for use in path regulation. 

In flight, these latter two controls are calculated by the trim solution 
algorithm to generate the total acceleration command, including the corrective 
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acceleration commands from the path error regulator, at the current configura- 
tion and flight condition. In order to maximize the available control margin, 
the trim solution algorithm also uses the nozzle for regulation if one of the 
principal regulator controls saturates. 


Report Outline 

The basic Trimmap design problem is formulated in the next two sections. 
In the second section, "Trajectory Commands," the input acceleration commands 
and trajectory variables are derived for the various maneuvers and steady 
flight conditions typical of flight operations; the domain for these input 
variables is estimated from the kinematic constraints imposed by passenger 
operations and by the limits of the AWJSRA performance. The AWJSRA force 
model and trim equations are given, including the domains of control variables 
and parameters for which the model is defined or which bound the region of 
design interest, in the third section. In the latter section, the nature of 
solutions of the trim equations at a given flight condition is also examined, 
including their redundancy and the effect of operational constraints on 
controls usage in limiting the range of redundant solutions which can be 
employed by the Trimmap. These constraints further define the operational 
flight envelope of the aircraft as those flight conditions for which solutions 
of the trim equations exist satisfying the operational constraints. From 
these considerations, the Trimmap design problem is formulated and a structure 
is outlined for the Trimmap which includes a configuration schedule to resolve 
the control redundancy and a trim solution algorithm for regulation. 

The trim solution algorithm is considered next. At a given flight condi- 
tion the acceleration capability of the aircraft for flight-path error regula- 
tion, or control margin, is found to be a function of the configuration and 
the controls used by ' the algorithm. A structure is outlined which maximizes 
the control margin with respect to the controls used by the algorithm. 

The configuration schedule, the main topic of this report, is developed 
in the next three sections. The schedule determines the operational flight 
envelope of the automatic system as well as the values of all configuration- 
dependent parameters of interest, such as control margin, attitude, fuel rate, 
etc., over that envelope. The development begins with a detailed examination 
of the operational constraints on controls usage for the AWJSRA and a selection 
of the parameters to be optimized by the schedule. Then, the flap schedule is 
calculated to provide maximum lift margin and the nozzle schedule is calculated 
to maximize control margin. This is done with an approximate numerical 
procedure which consists of enumerating solutions at a sufficient number of 
points covering the domains of interest of all the independent variables 
defining configuration and flight condition. The numerical results for the 
optimum configuration at each flight condition are fitted with approximating 
functions for use as the configuration schedule in simulation testing of the 
proposed control system. The analysis considers the flight envelope limits 
for the design and the behavior of the configuration schedule over that 
envelope and its effect on the remaining controls and other configuration- 
dependent parameters of interest. 
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Finally, the Trimmap is completed with the addition of elevator trim and 
attitude and angular rate command logic to the basic Trimmap and, together 
with the remaining control system of figure 1, is subjected to simulation 
testing on a STOL approach path. Behavior in response to various maneuver 
commands, Trimmap model errors, and wind disturbances is examined. 


TRAJECTORY COLIMANDS 


The nature of the acceleration commands to the Force Trimmap is examined 
in this section. Path axes are introduced as a convenient reference frame 
in which to define the acceleration commands associated with aircraft maneuvers 
and trajectories, and in which to state the operational constraints on trajec- 
tory kinematics associated with passenger operations. On the other hand, the 
basic Trimmap works with the stability axis components of the applied specific 
force. These are functions of the trajectory commands; their relation to the 
path axis components of acceleration, their behavior during typical aircraft 
maneuvers, and their admissible ranges of values for passenger operations, 
which forms the Trimmap design domain, are examined. 

The aircraft velocity vector, expressed in spherical coordinates, is 

V = V(cos Y cos i|j^i + cos y sin i[»yj - sin yk) (3) 

where (i, j, k) are unit vectors defining runway axes and (V, y) are the 

aircraft speed, heading angle, and flight-path angle (fig. 3(a)); they can be 
calculated from the runway axes velocity coordinates using: 


V = }/x^ + y^ + 
tan ^ ^ [-' JfjTf) 

sin y = - Y ^ [-iT/2,Tr/2] 

The acceleration vector for any trajectory can be derived from equation (3) as 

a = cos ym - Vyn (5) 

where m, n are unit vectors defining the path axes (fig. 3(b)) which are 
oriented, respectively, tangent to the flight path and normal to the flight 
path in the horizontal and vertical planes.^ They are related to runway axes 
by 


^The underbar, u, is used to distinguish the unit path tangent vector 
from the control vector, u, used elsewhere. 
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u = cos Y cos ilj^i + cos y sin - sin yk 
m = -sin i|;^i + cos ipyj 

n = sin Y cos i/;^i + sin y sin i/^yj + cos y^ 

The path axes components of a in equation (5) separate the maneuver 
acceleration command into terms corresponding to speed changes (V) , turning 
flight (V^/y cos y) or, equivalently, (V^/R^^) , and flare or pushover maneuvers 
to change flight-path angle (Vy) . Maneuvers from one steady flight condition 
to another can be formulated as a superposition of these maneuvers and, for 
CTOL aircraft, maneuver accelerations along the three axes are controlled 
independently by engine throttle, roll attitude, and pitch attitude. The path 
axes are unsuited to hovering (V = 0) and vertical flight (y = x = 0) , for 
which they are undefined. 

The aircraft trajectory is governed by 

ma = W + + Fg 

where the terms on the right-hand side are weight and aerodynamic and engine 
forces and m is the aircraft mass. The acceleration command of interest in 
the Trimmap design is the applied specific force to be supplied by the engine 
and aerodynamic forces that the Trimmap has at its control. Given in g-units, 
the applied specific force command is 




f = g) 

g 



+ V 


w 


(7a) 


_ •T' 

where g, referred to path axes, is g(sin y, 0, cos y) . The stability axis 
components of f are the most useful form for the Trimmap calculations and 
are calculated as the cylindrical coordinates of f (see footnote) : ^ 


A = u • f = sin y H 

u - ' g I 


= I [I - I = yj (cos Y + Vy/g)2 + v2/gR^ 

V2/gR, 


= tan 


\n • f / \cos y + Vy/g y 


(7b) 


^The derivation of equations (7) assumes zero nominal sideslip angle, 
so that the stability x-axis coincides with the velocity direction, u. While 
this is usually the desired value of sideslip can, in general, be 

commanded independently of the trajectory. General formulas, which include 
the effect of nonzero Sq used in the Force Trimmap, are derived in 

appendix D. 
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where is the longitudinal applied specific force, is the applied 

specific force normal to the flight path, and ())y is the angle of the normal 
specific force vector from the vertical plane. This angle is measured about 
the velocity vector and differs in definition from the usual Euler roll atti- 
tude angle measured about the body x-axis only by second-order effects. 

Equation (7) maps the trajectory command into corresponding applied 
specific force commands for the Trimmap and is used in figure 4 to map a 
graph of (A^, A]^) for the principal flight conditions of interest; for static 
equilibrium flight (a = 0) , equations (7) generate a unit circle as y is 
varied. For steady turns (|Rc 1 < ") > the unit circle is moved up and for 
steady speed changes (V ^ 0) these circles move left or right by the amount 
V/g. These constitute the flight conditions which can be held for a more or 
less extended period. Acceleration commands for a typical approach trajectory 
will consist of a sequence of such steady operating points together with the 
maneuvers required to pass from one steady operating point to another, such as 
the flare and pushover maneuvers shown in figure 5, and the example command 
generator output for part of an approach path shown in figure 6. 


The nominal path maneuvering commands to the Trimmap are limited by the 
aircraft capability, passenger comfort constraints, and operating constraints 
for the specific aircraft (ref. 3). For the AWJSRA these limits are taken as: 


h 

G 

[0,3500 

m] 

Va 

G 

[45,160 

knots] 

YaI 

< 

15° 


Ivl 

< 

0.15 g 


v;i 

< 

0.125 g 


Ul 

< 

20° 





> 


/ 


( 8 ) 


In the present context, these limits are assumed to be satisfied by the 
trajectory commands of both the ATC and the command generator. The limits on 
|Vy|, \(^\ typical of passenger operations, and the remaining are esti- 

mated outer bounds of the AWJSRA flight capability and correspond to wing 
stall and the extremes of total applied longitudinal specific force 

. From these and equations (7) the extreme nominal 
' 'min; 'max’ min’ max' ^ ^ ' 

specific force commands that need be considered in the present Trimmap design 
are therefore enclosed in a regime: 


e [-0.35,0.35] 

Aj, e [0.75,1.5] J 

as shown in figure 4. The most used portion of this regime for passenger 
STOL and CTOL operations is within ±0.2 g of level static equilibrium flight. 
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but other types of flight operations (helicopters, military aircraft, rockets) 
extend well beyond this region. 


AUGMENTOR WING TRIM EQUATIONS AND TRIMMAP FORMULATION 


Trim Equations 


The Trimmap uses the controllable engine and aerodynamic forces to gener- 
ate the commanded applied specific force. For STOL and CTOL aircraft these 
forces are given in stability axes and are assumed aligned in the aircraft 
plane of symmetry; for the AWJSRA these are expressed as 




/cos (a + v)v 




M . 

\si 


in(a + v), 



(9) 


The engine thrust is vectored through controllable exhaust nozzles (fig. 7) 
whose angle from the longitudinal body axis is v. In addition, compressed 
cold air from the engine is exhausted through a specially designed flap 
(fig. 7) to augment the aerodynamic forces. The potential thrust T^ of this 
compressed air, obtainable by isentropic expansion to the atmosphere, is used 
to define the cold thrust coefficient 


C 


J 



and the dependence of aerodynamic forces on engine power is given in terms 
of this coefficient; that is. 


C 


D 


C 


L 


— Cji^ (ot , 6^ , C j) 


} 


( 10 ) 


The engine parameters Tjj, Tq, mg are functions of the throttle setting, 
airspeed, and ratios of stagnation air temperature (x) ^and pressure (6) to 
reference sea level values: 


Tc = Tc(6t,V^,T,6) 


( 11 ) 


The engine and aerodynamic model parameters listed in equations (10) and (11) 
are nonlinear scalar functions available as tabulated data. The modeling data 
used for this study are given in appendix A, and are taken from a detailed 
Ames simulation of the AWJSRA based on wind-tunnel and flight-test data 
(refs. 4,5). 
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The most accurate available model data are desired as the starting point 
for the study. Model data are usually available from simulations; their 
accuracy will depend on the source (theory, wind-tunnel tests, engine test 
stand data, flight tests) and usually improves with time as more data becomes 
available. Errors in the static force model affect system performance as a 
result of inaccurate calculation of safety margins, optimal configurations, 
and system flight envelope; in flight, steady-state acceleration response 
errors are compensated by automatic adjustment of controls from their pre- 
dicted values and, if sufficiently large, induce steady position errors. On 
the other hand, increased complexity and computational requirements are 
associated with increased model accuracy and some simplification of the 
simulation model is desired both for off-line studies and for the on-board 
Trimmap. For the present study, the static force simulation model is repre- 
sented to 0.01 g by the model of appendix A. This value is small relative to 
the assumed compensation authority of the control system so that that authority 
will be used in flight on the residual uncertainties of the simulation model 
itself. 


Trim Equation Variables 

Equation (9) must be solved to obtain the output of the basic Trimmap. 
For this purpose it is useful to sort the variables of equation (9) into 

Control variables: u = (6£,v,a,6^) 

Path variables: x = (V£,A^,An) 

Parameters: p = (W,x,6) 

where Vg denotes the equivalent airspeed V^^(h)/p(o). The aircraft 
model in equation (9) is defined only over a bounded domain of values of u 
associated with hardware limits on the controls and the availability of model 
data. In the present case, the aircraft model is defined for 

= {u: e [5.6°,72°] , 

V e [6°,104°] , 
a e [-10.5°,27.5°] , 

6^ e [0°,min{38°,36.3 + 143. 5(/^ 

The variable a is controlled with the elevator, whose travel is limited to 

6e^ ^ [-25M5°] 

The required elevator setting to trim the aircraft at any value of a can 
be calculated from the pitching-moment model of the aircraft. Results have 
been obtained for this trim setting, dgrp (a , 6f ,Cj) , and are included in appen- 
dix A. 
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The path variables x are calculated from the commanded path. Their 
ranges of interest in the Trimmap are limited by aircraft capabilities, 
passenger comfort, and operational constraints as noted in the previous sec- 
tion. An outer bound on the regime of interest for the AWJSRA is estimated as 

^ [^5,160 knots], € [-0.35,0.35] ,Ajj e [0.75,1.5]} (14) 

Not all points in are necessarily within the aircraft capabilities 

defined by equations (9)-(13), and a more exact determination of these 
capabilities is part of the Trimmap derivation. 

The weight and atmospheric parameters, p, which appear in eq. (9) will 
also take on values over some desired domain, within which satisfactory 

Trimmap operation should be provided. For the AWJSRA this domain can be taken 
as 

^o='(p; w e [0.95,1. 2]Wg,T e [0.9, 1.1], 6 e [0.7, 1.0]} (i5) 

where Wg is the standard weight value 

Wg = 177.9 kN (40,000 Ibf) 

The ranges of the atmosphere parameters in ^ is smaller than for passenger 
operations generally and reflects the aircraft’s operation principally for 
research in a single terminal area. 

In summary, aircraft trim equations, such as equation (9), can be given 
in the general form 


a^(x,p,u) =0, U E 

where x, p, u refer, respectively, to path variables, parameters, and 
control variables. The operational domain over which its solutions are of 
interest in the Trimmap is bounded and given in the form 

P ^ ^ 

Because the trim equations are often nonlinear and lack closed-form solutions, 
numerical procedures are necessary in the Trimmap analysis and the magnitude 
of the task is sized by the dimension and extent of these domains. 


Control Redundancy 

The trim equations are to be solved for control settings u given the 
flight conditions (x, p) . There are, however, redundant solutions whenever 
the number of controls exceeds the number of equations. For the AWJSRA this 
excess is two. 


dim[u] = dim[a^] + 2 
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and there can be arbitrarily many solutions of equation (9) within the 
limits, The set of all such trim solutions at a given flight condition 

can be formally denoted as 

= {u e aY(x,p,u) = 0 J (16) 

and the aircraft flight envelope is then all those flight conditions within 
the design domain, for which there is at least one trim solution: 

= {(x,p): X e ^^,p e ^ and ^^(x,p) 0 J (17) 

Although there are no closed form solutions of equation (9) for ^«p(x,p) , 
numerical solutions can be generated in sufficient quantities to produce a 
plot of for any single flight condition. Figure 8 shows such plots 

for two flight conditions (level flight at 120 knots and descending flight at 
-7.5° glidepath angle and 65 knots). Every point in these figures yields a 
solution by interpolating the map lines for the flap and nozzle values at that 
point . 

Solution of the nonlinear trim equation will be required repeatedly in 
this analysis and in the actual flight Trimmap, so that efficient solution 
algorithms are significant to the derivation. Algorithms for solving the trim 
equation are described in appendix B. 


Operational Constraints on Controls 

Not all the trim solutions in ^rp(x,p) are acceptable for actual steady 
flight operations as Trimmap output commands; only those which satisfy a 
number of operational and safety constraints on the use of controls (limits 
on throttle, lift margin, angle of attack, etc.) are acceptable. This reduced 
set of acceptable trim solutions is denoted by 

^ ^ {fi(u>x,p) >0, 1 = 1,..., N}} (18) 

where (fi) are all the operational constraints that the Trimmap output must 
satisfy in flight for the particular aircraft. These are listed in detail for 
the AWJSRA in the next section. The effect of imposing these constraints 
on the trim solution sets previously given for two flight conditions is 

shown in figure 9. As seen, is considerably reduced from in both 

cases. 

Operational criteria for powered-lift aircraft are less well defined 
than for CTOL aircraft and are the subject of recent and ongoing studies at 
Ames (refs. 6-8). The criteria used in the present study are taken from 
current AWJSRA research flight operations and, in their details, are unique 
to the AWJSRA. However, constraints of the type given and their effect on 
the design method would occur for aircraft generally. A minimum requirement 
for the constraint list used in the computations is that it comprises suffi- 
cient conditions for operational acceptability of trim solutions. 
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The operational flight envelope for the aircraft within the design domain 
is that set of flight conditions for which at least one trim solution accept- 
able for steady flight is available; that is, 

= |(x,p): and o} (19) 

Thus, the constraints define the limits of operationally acceptable flight as 
well as restrict the acceptable trim solutions at any point within 


Trimmap Formulation 

A functional description of the Trimmap element of figure 1 can be given 
in greater detail. The Trimmap function is to solve the aircraft trim equa- 
tions for a unique control coimand, u^, given the path commands and 

current parameter values p. Over the operationally acceptable flight envelope 
jr,pA it resolves the control redundancy by selecting an optimum trim solution 
subject to constraints on controls usage. If the Trimmap obtained in a design 
is denoted as the map 


u*(x,p): ^ ^ (20) 

then that design will achieve some corresponding operationally acceptable 
flight envelope 

•^TA ^ ^ ^ -^’P ^^>u*(x,p) ^ ^^^(x,p)} (21) 

and complete flight envelope 

E {(x,p): X ^ ^ ^,u*(x,p) e^j^(x,p)} (22) 

At best, these are equal, respectively, to and ^ (eqs. (17) and (19)) 

but are usually smaller than those two sets as a result of approximations and 
simplifications made in the design process. An objective of the Trimmap 
design is to maximize and since these are outer limits on the flight 

capabilities of the autopilot-aircraft combination. 

The basic structure of the proposed Trimmap logic is shoxm in figure 10. 
Two principal elements are used: a configuration schedule and a trim solution 

algorithm. This structure assumes that an optimum trim solution satisfying 
the operational constraints for steady flight, u*(x,p), has been calculated 
off-line and solutions from u* for two of the controls, (u^, u^) , are stored 
on-board as the configuration schedule. In flight, the configuration schedule 
commands these two controls on the basis of the nominal trajectory commanded 
by the command generator; that is, these two controls are driven open-loop by 
the nominal trajectory command. The remaining two controls, (u^, u^), are 
calculated by solving the two trim equations using the current parameter 
values p and x^, where x^ differs from Xq by the amount of the corrective 
acceleration command from the trajectory regulator. These two controls are in 
the feedback line and are termed regulator controls here. Their command values 
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will differ from their nominal values in u*(xo>p) to the extent that path 
errors are nonzero. 

Further details of the above structure applied to the AWJSRA are shown 
in figure 10(b). Flap and nozzle have been selected as configuration controls. 
Their stored schedules, F*(x) and v*(x) do not include the parameters p as 
independent variables but, instead, the configuration schedule is optimized 
for a standard set of values, Pg, given by 

Wg = 177.9 kN 



This simplification is reasonable in view of the narrow regime of parameter 
values ^ over which the AWJSRA is operated. In flight, the trim solution 
algorithm adjusts the regulator controls for nonstandard parameter values, but 
some reduction in the design envelope is traded for the simplification. 

Of the four controls, only the flap actuator response is too slow for use in 
path regulation but any pair of controls from {v, ot, 6^-} or even all three 
controls can be used as regulator controls. The pair (a, 6^) is used in CTOL 
aircraft and was also found suitable for AWJSRA operations as the basic 
regulator mode throughout its flight envelope. However, the use of all three 
controls provides greater regulator control power than is available for any 
single pair, and the trim solution algorithm is modified to utilize the nozzle 
control when the basic regulator mode saturates, as discussed in the next 
section. 

The basic Trimmap outlined in figure 10(b) is that element of the Force 
Trlmmap of figure 2 that deals with the solution of the two trim equations, 
equation (9). This element encompasses the principal design problems regarding 
nonlinearity and trim solution redundancy. Additional functions are required 
in the Force Trimmap to map input acceleration commands into stability axis 
components (eqs. (6) and (7)), to provide the elevator trim setting that 
balances out pitching moments at the commanded configuration (appendix A) , and 
to provide appropriate attitude commands that orient the applied forces 
correctly in inertial space (appendix D) . 


STRUCTURE OF THE TRIM SOLUTIpN ALGORITHM 


Three of the AWJSRA’ s controls are suitable for path error regulation. 

In general, the control margin available for the regulation function differs 
with the control pair used and the three controls (6t, v, 6^) collectively 
provide more control margin than is available from any single pair. Thus, 
to obtain maximum control margin the trim solution algorithm must use all three 
controls. This section includes a discussion of the regulator control margin 
and outlines a regulator control logic for the AWJSRA that uses all three 
controls . 
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The acceleration capability or control regime of the regulator controls 
at a given configuration {6f, v} and parameter values p can be calculated 
as those accelerations that can be generated in equation (9) by moving the 
regulator controls over their permitted ranges of use for regulation; that is, 
if {a, 6|-} are the regulator controls (fixed nozzle mode) then the control 
regime is the envelope of accelerations filled out by 


(23) 

It is understood in equation (23) that {<Sf, v, Vg, p} are given and that 
{ct, 5^} are varied through their permitted ranges for regulation, 

• The numerical values of these limits on control 

t^in 4nax ’ 

usage are lisited in the next section along with other constraints. The con- 
trol regime is illustrated in figure 11 for glide slope conditions. 

Control regimes for other regulator control pairs or for all three 
controls can be defined analogously 


" {(Au.An); V £ ^ 


mm 


^5^(Ve,P) = {(A^.An); V e [Vj^in.^max] »<5t ^ ’ 

a ^ [a„. ,i„(5c,p,u) = O} 

min max i > 


(24) 


These different control regimes are compared in figure 12 for glide slope 
conditions; it is apparent that the control regime for all three controls 
includes and is significantly larger than the regimes of any of the two-control 
regulators. 

A method of measuring the control margin associated with such irregular 
figures is necessary for the configuration optimization analysis of this 
report which attempts to maximize the control margin at all operating points. 
The control margin at a nominal operating point (A^^, A^q) is defined as the 
radius of the largest ellipse of the form 


r2 = ct2(^ _ _ Aj^^)2 (25) 

which can be inscribed in the acceleration envelope. This is illustrated by 
the ellipse in figure 11 which is centered at the glide slope operating point. 
The definition of control margin contains a selectable parameter o which is 
the axis ratio of the ellipse. It represents a weighting of normal accelera- 
tion relative to longitudinal acceleration in measuring control margin and 
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was selected for this study as the ratio of minimum acceptable acceleration 
margins for the normal and longitudinal axes 

0 = 5 

This value reflects the approximate ratio of disturbance accelerations and 
corresponding corrective acceleration activity of the path regulator for the 
two axes which results from the high lift-drag ratios, on the order of 
5 to 10, for aircraft. To formalize the definition of control margin, note 
that the region occupied by an ellipse of the form of equation (25) is 


and the control margin at a given flight condition and configuration is then 






(2,7) 


A method for generating the figures etc., and calculating control 

margins is described in appendix C. 


Control margins for the other regulator modes and for the three-controller 
mode are defined analogously 

\ 


(28) 



= maxj 

[R: 


(^>^Uq »^No) 

C 

•^6 


= maxj 

[R: 



C 



= max- 

{R: 

S 1 


C 

•^0 


Values of these control margins, measured in g*s, are noted in figure 12 for 
the various modes and it is apparent that the three-controller regulator has 
the largest margin. 


A conclusion drawn here is that the regulator logic must use all three 
controls in order to realize the maximum available control margin. One method 
of implementing this for the AWJSRA is to use the fixed nozzle regulator near 
the nominal operating point but to switch to an alternative mode which varies 
the nozzle whenever one of the controls {a, is saturated at an extreme 

permitted value. A corresponding mode-switching diagram for the glide slope 
case is shown in figure 13; the control regime subdivided into 

regions in which one or another of the three controls is held fixed. A flow 
diagram for the corresponding trim solution algorithm is shown in figure 14. 


The fixed-nozzle mode is selected in the present study as the control 
mode used near the nominal operating point, or central regulator mode. The 
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configuration optimization will be aimeH in part at maximizing the control 
margin of the central mode in order to minimize control activity and the 
likelihood of saturating the controls of the central mode. The fixed nozzle 
mode uses the same pair of control variables as is used for regulation in CTOL 
aircraft and is obviously appropriate for the AWJSRA. in its conventional 
cruise configuration. The present Trimmap extends this mode throughout the 
AWJSRA flight envelope and attempts to preserve good control power and accept- 
able sensitivity into the low-speed STOL regime and powered-lift configuration. 
It does this by selecting the configuration to maximize the control margin of 
the central mode and by resorting to the use of the third control when the 
central mode saturates. More generally, alternative choices of central mode 
can be considered in regions of the flight envelope where improved control 
margin is both needed and available from an alternative mode. It should also 
be noted that, while the multiple mode structure of figure 14 realizes the 
maximum available control margin, further improvements in control sensitivity 
should be available with an algorithm which uses all three controls simulta- 
neously for regulation, but such an algorithm is outside the scope of the 
present work. 


CONFIGURATION OPTIMIZATION 


The configuration schedule is to be derived off-line and stored on-board 
and the derivation seeks to optimize the configuration subject to the opera- 
tional constraints on controls usage for the aircraft. This section describes 
the AWJSRA constraints and the optimization procedure used. 


Operational Constraints for the AWJSRA 

The operational constraints which bound the acceptable control settings 
for steady flight and for regulator usage are discussed below and listed in 
figure 15(a). These are limits on the extremes of each control or on functions 
of the controls. Many details of the constraints are unique to the AWJSRA 
and its operations as a research aircraft, but the list is expected to be 
representative of most powered-lift aircraft in terms of length, complexity, 
and type of functions to be bounded. Criteria for civil operation of powered- 
lift STOL aircraft are under development (refs. 6-8). 

The flap range is bounded by hardware limits and by a structural placard 
on the maximum value which is a function of equivalent airspeed (fig. 15(b)): 
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6f (Vr) = mini65,max{3(106.7 - Ve) ,0.66(165.5 - Vr),5.6}J (29)'* 

max ' ' 


Throttle servo stops have been installed on the aircraft to limit the 
extremes used during automatic flight control; the minimum stop is operated 
by the flap setting and is satisfied in the present design by the continuous 
piecewise linear function (fig. 15(c)): 

= min{l7.6,max{0.553(6f - 13. 2), 9. 3}} (30) 


This minimum is well above the engine idle setting and is selected to limit 
the emergency power recovery time from minimum throttle. Shorter recovery 
time is required for flaps-down powered-lift configurations than otherwise. 

The maximum throttle setting for nominal or steady flight conditions is 
equivalent to the maximum continuous power setting established for the engine: 


^max , NOM 


25.6° 


(31) 


This is the upper limit imposed on throttle settings in the configuration 
schedule derivation. Temporarily higher settings are permitted for the regu- 
lation function and for takeoff; this upper limit is implemented by one of the 
servo stops. In addition, the maximum throttle is further reduced on cold 
days due to a structural limit on the thrust force which the engine nozzles 
can support when nozzles are down. Combining these limits gives (fig. 15(d)): 

•Sf = mini 29. 4, 136. 6/r (0.705 + 0.295/6) - 107.174} (32) 

unax j REG 

Nozzles are limited simply by their hardware extremes as noted in figure 15(a). 


The lift margin at a given flight condition (x, p) and control setting u 
is the amount of lift force acceleration that can be generated by increasing 
a to its stall or maximum permitted value; that is, 


LM(x,p,u) E [Ci^3^(6f,Cj)-CL(a,6f,Cj)] 


(33) 


The lift margin reflects both a-margin from stall and speed-margin from its 
minimum value for equilibrium flight at the given flap, power and nozzle 
settings; therefore, it measures the tolerance of the aircraft at (x, p, u) 
to vertical gusts, windshears and speed control errors as well as the ability 


^The use of minimum and maximum functions, as in equation (29), provides 
a conveniently compact notation for representing piecewise linear functions in 
one or more variables. Functions of one variable, such as equation (29) are, 
however, more easily visualized in the standard notation; for example. 


max 


(Vr) 


f65 

3(106.7 - Ve) 
0.66(165.5 - Ve) 
< 5.6 


Ve < 85 
Ve e [85,90] 
Ve e [90,157] 
Ve > 157 
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to flare the aircraft from (x, p, u) using attitude changes only. A minimum 
value of lift margin is required at all nominal flight conditions and, for the 
AWJSRA, suitable margins established from flight experience are 0.69 g in the 
flaps-up configuration and 0.4 g in flaps-down (6f > 30°) powered-lift config- 
urations. These requirements are satisfied by the constraint: 

LM(x,p,u) > max *10.4, min{0. 69, 0.0119(63. 66 - 6f ) }| (34) 

The value of 0.69 g matches the usual speed margin of 30% required for CTOL 
aircraft. In a powered-lift configuration, however, lift can be generated 
from engine power increases as well as from angle of attack increases and a 
lower fixed-power lift margin can be comfortably permitted provided a fixed- 
power flare is possible (ref. 7). 

The lift-margin constraint limits the acceptable control settings at a 
given flight condition. It can also be viewed as placing an upper bound on 
the nominal value of a as a function of {5^, 6^, x, p} or as establishing the 
the minimum operating speed boundary of the flight envelope. This boundary 
is shown in figure 15(e) for values of (6^, 6^, x, p} representative of the 
approach flight regime. 

The constraints on pitch attitude listed in figure 15(a) are approximate 
limits favored by pilots. They are imposed in computing the configuration 
schedule by corresponding approximate constraints on a(x, p) ; assuming 

1 . a = 0 - Y 

2. Y = sin“hA^) 

Then the limits on a for nominal flight conditions are 

“NOM ^ [-10'' - Sin"^ (Ay) ,15° - sin"^ Au] 

I 

Secondary effects of roll angle are neglected in (1) above and y has 
been calculated from A^ assuming V is null. More generally, if V 
nonzero then, in the extreme cases, equation (35a) implies 


1. 

If 

a = 

“NOM.min 

then 0 > 

for V < 0 


2. 

If 

a = 

°^N0M,max 

then 0 < 0 jjj 3 ^ 

for V ^ 0 


and this ensures 

the 

pitch attitude constraints 

in most operational 

situations 

with nonzero 
by requiring 

V. 

In 

addition. 

the constraint on 

touchdown attitude 

is imposed 


(35a) 

is 


°^N0M ^ ^ sin ^ (A^) (35b) 

at those values of {A^, Vg} for which touchdown can occur. 

In flight, the regulator varies angle of attack around its nominal value; 
an upper bound to prevent stall by the regulator is established as 
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“max, REG " min|0.1667 (Vg + 10), 15} 


(36) 


This limit is used in calculating the control margin of a given configuration 
and provides a margin from the stall angle of attack that gives at least 
0.15 g of lift margin unused by the autopilot in the case of approaches at 
maximum flap setting. These boundaries are shown in figure 15(e). 


In flight, path errors are corrected by varying the regulator controls 
around their nominal values given from the configuration schedule. Conse- 
quently, the choice of configuration controls must allow sufficient margins 
from the boundaries of aircraft acceleration capability for use by the 
regulator controls. Lower limits on these margins were taken as 


Longitudinal acceleration margin ±0.05 g 

Normal acceleration margin ±0.25 g 


(37) 


From equation (7) and the lift force expression it can be calculated that 
these values allow for corrective flight-path variations of ±3® or speed 
corrections at ±0.05 g and can accommodate horizontal gusts of 13%-15% of the 
airspeed without exceeding the margins. In terms of the control margin measure 
defined in equations (25) and (28), the above acceleration margin limits are 
given by the constraint 


CM5^(xq,p) > 0.25 g 

This is the minimum acceptable margin when all three regulator controls are 
considered. The control margin of the central regulator mode, which uses only 
two controls, can be less than this value. In the glide slope case previously 
shown in figure 1,2 the desired control margin is available using all three 
controls but not for any choice of two-controller modes. 

The elevator is used to execute angle of attack commands from the Trimmap 
and to regulate pitch attitude errors. A trim elevator setting associated with 
the nominal flight condition, and calculated as that elevator setting that 
trims out all pitching moments due to engine and aerodynamic forces, is derived 
in appendix A. In the present configuration schedule computations a margin 
of 8° from the limits of elevator travel was imposed on the associated trim 
elevator setting 


6e^(x,u) G [-17°, 7°] 

in order to provide some margin for the regulation functions. 

The configuration schedule is derived by identifying the operationally 
acceptable range of choice of configurations for steady flight (eq. (18)) 
and then determining a unique choice that optimizes some parameter of interest. 
The range of acceptable configurations is defined by constraints on controls 
usage, such as those listed in figure 15(a) for the AWJSRA. These constraints 
considerably reduce the range of choice compared to the full set of trim solu- 
tions and also impose gross variations of configuration over the flight 
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envelope. Both of these effects can be seen in the two examples previously 
illustrated in figure 9. 


Configuration Selection 

Schedules for the flap and nozzle settings are to be determined over the 
flight envelope. Schedules that optimize control characteristics are computed 
next in two successive steps using the criteria: 

1. Select the maximum acceptable flap angle 

2. Select nozzle angle to maximize the regulator control margin 

The maximum flap is selected within the acceptable range of choice and 
can be calculated independent of any particular choice of nozzle schedule. 

This flap selection criterion maximizes the lift margin, which is a strong 
function of flap. It emphasizes safety and control in the low speed approach 
regime and is also favorable for fuel requirements in this regime. This occurs 
because flap and throttle provide lift margin interchangeably there so that a 
higher flap setting permits a lower range of throttle settings for trim 
(e.g. , see fig. 9(b)) . 

After the flap schedule is determined, the nozzle is selected to maximize 
the control margin of the central regulator. This margin depends strongly on 
the nozzle setting at any given flight condition and flap setting. This 
selection criterion minimizes control sensitivity and saturation of the central 
mode, but does not affect the control margin available for all three controls, 
CM^^. 


An additional step is useful to adjust the schedules derived from the 
above procedure. After the variations of all parameters over the flight 
envelope 


F’^(x), v^(x), 6*(x), a*(x), . . . 

are obtained, the scheduled controls can be adjusted for any engineering 
requirements for smooth or monotonic variations, or for special restrictions 
on controls in restricted regions or for certain maneuvers. 

The configuration selection procedure is carried out empirically in the 
next two sections; the magnitude of the computational task is reduced in the 
stated procedure by optimizing the flap and nozzle schedules in sequence 
rather than simultaneously; this is based on the special properties that the 
parameters to be optimized are each strongly dependent on only one of the con- 
figuration controls, and that the flap schedule can be optimized independent 
of the nozzle setting. These specializations reduce the domain of search in 
the optimization; but their effect on the margins obtained is expected to be 
small. 

The configuration selection rules used emphasize control characteristics 
on the assumption that lift margin and control power are marginal in much of 


22 



the low-speed STOL approach regime. A major objective of the configuration 
scheduling is then to maintain acceptable control characteristics and to 
realize as much of the potential operating envelope of the aircraft as is pos- 
sible in this marginal regime. A second generation design can be expected to 
integrate fuel and noise considerations into the optimization by trading excess 
lift or control margin obtained in the present design in some regions of the 
flight envelope for reduced fuel rate and noise. 


FLAP SCHEDULE 


The flap schedule to be used in the ACS is described in this section. 
The computations are done in two steps; first, the range of admissible flap 
settings is determined at each of a number of operating points covering the 
trajectory envelope of interest, and then a unique setting is defined at 

each point and a function fitted to these values for use in the ACS design. 
The flap schedule given here is approximately the maximum admissible flap at 
each point in The influence of flap servo rate limits on the design 

of the flap schedule is also analyzed. 


Admissible Flap Settings 

At a given flight condition, (x, p) , the admissable flap settings are 
those values for which at least one acceptable trim solution exists. An 
algorithm that was used to calculate these flap settings at a sufficient num- 
ber of points is outlined in appendix C. The procedure is entirely empirical; 
a grid of values of flap settings is tested at a given flight condition and 
then this is repeated for a grid of values of (x, p) covering the flight 
conditions of interest. 

Results are shown in figures 16(a)-(g) as a series of plots of the admis- 
sible flap range for several positive values of A^ (climbing or accelerating 
flight) and negative values of A^ (descending or decelerating flight) . The 
standard parameter values and a value of 1 g for A^ were used in these 
calculations. The upper bound on the admissible flap range is limited by the 
flap placard and, at high values of A^, by the control margin. The lower 
bound is dictated at low speeds by the lift margin constraint. 


Flap Selection 

The rationale for flap scheduling is given below and is basically a 
generalization of conventional flap usage. First, the flap schedule will be 
given as a function of only two of the variables of (x, p) : 

F* = F*(Ve,A^) (38) 

This choice provides for the major trends in conventional flap usage; that is, 
flaps are lowered with decreasing speed in order to maintain lift margin with 
the higher maximum lift coefficients of the flaps-down configuration. Second, 
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flaps are raised with increasing flight-path angle (or A^) and with increasing 
speed in order to avoid the induced drag penalty of the flaps-doxm configura- 
tion which must be overcome by engine thrust. 

Third, the flap will be scheduled to give the maximum admissible flap set- 
ting, or equivalently, to maximize the lift margin. This choice matches cur- 
rent practice with the AWJSRA in the approach regime and minimizes angle of 
attack with respect to flap setting. It is unfavorable for fuel requirements 
at higher speeds but favorable at lower speeds so that a minimum fuel design 
would differ from the present design in some regions of the flight envelope. 
Lift margin is an appropriate control characteristic to optimize in the flap 
selection since it is a strong function of flap. Optimization of the control 
margin, CM^^, was also considered since it is a function of the flap angle 
and is independent of the nozzle selection, but it was found that CM 5 £ was 
insensitive to flap setting. 

The effect of flap choice on the remaining controls of the trim solution 
is illustrated for two flight conditions in figure 9. The flap selection 
removes one degree of redundancy and narrows the remaining choice to a single 
constant flap line in figure 9. In the glide slope case the maximum flap 
provides the widest choice of acceptable trim solutions, the minimum throttle 
and fuel rate, and the minimum ot and maximum lift margin compared to other 
flap values. In the higher speed case the lowest throttle setting occurs at 
minimum flap. In both cases the usual interchange between flap and angle of 
attack as lift generators is evidenced by the nearly parallel flap contours 
and, more generally, is present independent of configuration and flight condi- 
tion. At low speeds engine power is also used to generate lift and a similar 
flap-throttle interchange occurs while at high speeds the conventional relation 
between flap and throttle setting due to the induced drag variation with flap 
is obtained. 


Flap Schedule 


The flap schedule will be given in the form of a piecewise linear function 
suitable for use in the Trimmap. This function was derived to fit the maximum 
admissible flap within the domain of x for which a nonzero range of admis- 
sible flap settings is available; outside this domain the fitted function 
extends linearly. The errors in fitting the maximum admissible flap are 
designed to maintain the scheduled flap inside its admissible range over as 
much of the flight regime as possible. The schedule is 


where 


F*(Au,Ve) 


Fi(Au.Ve) 

Fgo(Au) 


= min|FuB(Au) .max{F^ (A^,V e) ,5.6}} 


= min{65,max{538. 5 (0.1207 - Au),30}} 

{ 140 - Ve ^ knots 

Fgp(A^) + 0.6627(90 - Ve) Ve > 90 knots 

= min)50,max{333.3(0.215 - Au),30}} 


> ( 39 ) 
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In form, F* is a generalization of the flap placard function (eq. (29)), 
which is a linear function of Vg limited by the flap extremes; in F*, the 
parameters of that form have been made functions of 

The variation of F* with Vg is given in figure 16(h) for several 
breakpoint values of A^; the parameters defining this variation are linear 
functions of ^u between these values. Plots of F* are also superposed on 
the admissible flap ranges in figure 16 (a)- (g) and it is seen that there is 
almost no loss of acceptable flight envelope by the fitted function, F*. The 
parameters of F* are based principally on the standard case data (fig. 16) 
and therefore fit that case best. However, they have also been adjusted to 
maximize the acceptable flight envelope for values of A^ other than 1.0 and 
in order to adjust the implied pitch attitudes in some regions of the envelope. 

The complete flap schedule is shown by a contour plot in figure 17. The 
expected trends are visible; that is, F* decreases monotonically with increas- 
ing Vg or A^. The variation with Vg is independent of A^ for negative 
values of A^. Thus, pull-up to level flight from any steady descent is 
executed without flap configuration changes at any speed. For positive values 
of A^ flap reduces with A^ to a speed-dependent minimum value. The flap 
schedule F* is defined at all flight conditions both inside and outside the 
acceptable flight envelope, as seen in figure 17, even though it is fitted 
only to data inside that envelope. Inadvertent excursions of the aircraft 
outside the acceptable flight envelope enter a marginal region in which trim 
solutions exist but are unacceptable for steady flight. Owing to its linear- 
ity, the schedule, (eq. (39)), provides a marginal region of adequate size to 
accommodate such excursions. 


Flap Command Rates 

In the design process discussed in the previous section, and in the design 
of F* so far described, the scheduling problem is treated as one that can be 
solved independently at each point in the fligljit envelope assuming that each 
point is a steady or quasi-steady flight condition. The transient variations 
of Uc necessary to go from one steady flight condition to another were not 
considered nor were control rate constraints imposed on the schedule. However, 
the control servo actuators are usually rate or bandwidth limited and transient 
control lags will occur during maneuvers if this limit is saturated. In 
general, if x^(t) is the commanded history of flight conditions for any 
maneuver or trajectory within the domain of interest, and the Trlmmap is 
represented as u*(x), then the corresponding commanded control rates are 
given by 

Uc = 

These rates depend on both the maneuver and the Trimmap’s configuration 
schedule. 

The AWJSRA flap is the slowest of the controls and has the rate limit 

|6f| = 3.5°/sec (40) 

‘ ^ ' max 
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The flap command rates in maneuvering flight implied by the flap schedule are 


dt 


SF* . ^ 9F* • 
V + A, 


3Ve 


9Au 


(41) 


The partials of F* are piecewise constant and are shown in figure 18. The 
maximum maneuver command rates from the command generator are obtained from 
equations (7) and (8) as 

|v| < 0.1 g 


^ lYl 




0.125 g 
V 


For speed-change maneuvers the flap rate limit is saturated if 


96f^ 


f ' max 


V 


= 1.8 deg/knot 


max 


Since |3F*/3 Ve| is less than this value at all flight conditions (fig. 18(a)), 
the use of F* will not saturate the flap rate limit as a result of speed 
changes. The rate of change of is due primarily to commanded flight-path 

angle rates. Since y is confined to small angles, maximum values of A^ 
can occur only briefly and some flap rate saturation can be tolerated during 
flight-path angle changes. The value of the partial for which rate saturation 
occurs is 


3 6f 
^c 

3Au 


t ‘max 


1^1 


max 


Ve deg/g 


The partial of F* exceeds this limit by factors of 3 to 10 in those regions 
where the partial is nonzero (fig. 18(b)). Flap lags as high as 15° have 
occurred in simulated maneuvers. However, the Trimmap is designed to avoid 
maneuver execution errors during configuration lags by using the actual flap 
value in the trim solution algorithm rather than the commanded flap (see 
fig. 10(b)). This provides automatic adjustment of the regulator controls to 
compensate for flap lag so that the maneuver acceleration commands are 
accurately followed. 


To summarize, configuration control rate limits are of consequence only 
during maneuvering, 

Xc / 0 


since this is the only occasion for significant nominal control activity, 
u^ ^ 0. The effect of rate-limit saturation is control lag and a possible 
loss of acceptability of the resulting nonscheduled configuration. This 
should be prevented for long duration maneuvers by limiting the appropriate 
schedule gradient in accordance with maneuver limits of the command generator 
and the control rate limit. For short duration maneuvers configuration lag 
can be temporarily tolerated but must be correctly compensated by the regula- 
tor controls in order to obtain adequate tracking of the maneuver commands. 
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NOZZLE SCHEDULE 


The nozzle schedule to be used in the Trimmap is derived in this section. 
The nozzle is scheduled to maximize the control margin of the fixed-nozzle 
regulator, CM^^ everywhere in the acceptable flight envelope. A computer- 
generated solution for the optimum nozzle, v*(x), is given first and then a 
fitted function to be implemented in the airborne computer, v*(x), is given 
along with contour plots of the corresponding values of the other controls and 
parameters of interest which depend on the configuration schedule. The 
optimization calculations are carried out at the standard parameter values. 


Optimum Nozzle Angle, v* 

The computational procedure used to generate the nozzle setting for maxi- 
mum control margin is outlined in appendix C. The effect of nozzle setting on 
control margin is illustrated in figure 19 for the glide slope case; the 
control regimes for several values of nozzle angle along with the glide slope 
operating point are shown. The control margin at a given nozzle angle is 
related to the minimum distance from the operating point to the boundary of 
the corresponding control regime; as seen in figure 19, this distance is 
strongly dependent on nozzle angle, and is even undefined at the highest and 
lowest nozzle angles for which the control regimes do not enclose the operating 
point. 

The effect of nozzle angle on control margin is seen more generally in 
figure 20 for various values of A^ at low speed and at cruise speed (65 
and 160 knots) . The curves are dashed outside the range of acceptable trim 
solutions which is limited by bounds on lift margin, throttle, or attitude. 

In the low speed case (fig. 20(a)), definite control power maxima occur at each 
A^, and the optimum nozzle angle increases with decreasing A-u • The conven- 
tional minimum-nozzle configuration is not within the acceptable range at any 
value of A^ at this speed so that the aircraft requires vectored thrust to 
fly acceptably at this speed. In some cases the maximum occurs outside the 
acceptable range due to insufficient lift margin and a constrained optimum 
nozzle value with lower control margin must be selected. The corresponding 
throttle is seen to be the minimum for the acceptable range of solutions and 
the optimum solution at this speed has the minimum lift margin. 

At cruise speed (fig. 20(b)) there is considerably greater control margin 
than at 65 knots, but the curves again have definite maxima. These occur on 
the nozzles-up boundary at level and climbing flight and at higher nozzle angle 
for all descending or decelerating flight. The unconstrained optimum nozzle 
is now within the acceptable range and can be selected. However, the corre- 
sponding throttle setting is higher than at lower acceptable choices of the 
nozzle angle and, at this speed, some control margin could be traded for lower 
fuel rates by selecting lower nozzle angles. The required throttle increases 
with nozzle angle at all flight conditions shown in figure 20 and this relation- 
ship between throttle and nozzle occurs generally at other flight conditions. 
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The use of thrust vectoring is not actually required at cruise speed by the 
operational constraints except for very steep descent or to generate maximum 
deceleration. 

A contour plot of the optimized nozzle angle v* is shown in fig- 
ure 21(a). This is plotted for normal specific force at 1 g, which corresponds 
closely to static equilibrium flight conditions. The contours are roughly 
parallel and vary principally with longitudinal specific force. The nozzles- 
up conventional configuration is used for level and climbing or accelerating 
flight at speeds above 90 knots. Below 90 knots and in descending or deceler- 
ating fligjit the nozzles are lowered for maximum control power and reach near- 
maximum values at the extreme decelerating boundary. The boundary of this 
plot encloses the acceptable flight envelope within which v* is defined. 

In some regions of this envelope, v* is given by the unconstrained maximum 
control margin but otherwise is at a boundary of the range of acceptable nozzle 
angles so that some constraint parameter is at an extreme permitted value. A 
subdivision of this acceptable flight envelope into such regions is shown in 
figure 21(b). The critical low-speed region is characterized by minimum lift 
margin, while the minimum flight path angle region is characterized by minimum 
pitch attitudes. At the edges of the acceptable flight envelope a second 
constraint parameter, usually the throttle, reaches an extreme and the 
acceptable range of nozzle angles narrows to a single value. 

Fuel minimization is not considered in the present maximum control margin 
design of the nozzle schedule, but the fuel requirements of the schedule can 
be evaluated. Fuel requirements vary with nozzle setting since the required 
throttle setting and, therefore, fuel rate and noise, vary directly with 
nozzle angle (fig. 20). The acceptable nozzle settings for maximum control 
margin and for minimum throttle setting coincide in the two regions of con- 
strained optimization corresponding to minimum nozzle and minimum lift margin, 
but not in the region of unconstrained optimization. Consequently, some fuel 
reduction is possible in this region by trading excess control margin for 
reduced throttle settings. Data for such a schedule were computed and the 
computations showed fuel reductions over a regime of higher speeds and descend- 
ing flight (Vg > 110 knots, Au < 0) with maximum reductions in throttle setting 
at steep descents (y -- -5° to -7®). However, only the margins of this region 
are used by typical approaches so the nozzle schedule for the present tests 
was not modified to incorporate these improvements. 


Nozzle Schedule 

The optimized nozzle schedule is to be implemented in an airborne digital 
computer. Several methods are available, the choice being influenced by con- 
straints on storage and computation time in the airborne computer and by the 
ease of adjusting the implemented schedule in different regions of the flight 
envelope to iteratively improve the fit, or as may be required following 
simulation or flight tests. The available methods for modeling such empiri- 
cally determined functions include (1) programming the generating algorithm 
on-board and generating v* on-line, (2) storing the tabulated data 
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|'^IJk| ^ • • • 9 Nj, j 1, . . ., Nj, k 1, . . ., Nj^l 

that would be generated off-line and use table interpolation on-board, and 
(3) fitting an approximating function, v*(x), to the given tabulated data. 

Method (3) will be used here; method (1) cannot be used for lack of 
both storage and computation time in the airborne computer and method (2) 
requires large storage, in the amount Nj • Nj • but is easily adjusted. 

Method (3) requires the storage of relatively few parameters but the fitted 
function v*(x) is difficult to generate or adjust. 

The fitting function was selected and the parameters of the fit evaluated 
"by hand" using the procedure outlined next. The procedure is organized to fit 
the variation of with its independent variables, {A^, Vg}, one 

variable at a time and in the order of their importance in influencing the 
observed behavior of v* ; that is: 

1. Select a functional form, v*(c, A^) , to fit the observed variation of 
V* with A^. Here, c refers to the parameters of the function. For example, 
a bounded straight line fit is 

V* = + bA^,v„,in]] 

2. Select Vgj, Ajqj^ from {Vgj}, {A]^^} and calculate the parameters, 

cj^, which best fit the data i = 1, . . ., Nj}. Repeat this step for 

all velocities {Vgj} to obtain the set of parameters j = 1, . . ., Nj}. 

_3. Select a functional form, c(dj^, Vg) , to fit the observed variation 
of c with Vg. Here, d^^ refers to the parameters of the function. 

4. Calculate the parameter values, d^^, which best fit the data 

{cj^, j = 1, - . Nj}. Repeat these computations, beginning with step 2, 
for all values of normal specific force, to obtain the set of parameters 

{d|^, k — 1, • • ^ • 

5. Select and fit an approximating function, d(A^), to these data. 

This procedure res^a^ts in a function of the form 

v*(x) = v*(A^,c) 

where 

c = c(d,VE) 

5 = d(AN) 

The nonunique choices of fitting function and parameter values at each 
step of the above procedure allow useful latitude to impose engineering 
judgment. However, it has proved cumbersome, and a more algorithmic approach 
would be necessary if more than three independent variables were desired in 
the schedule. The functional forms selected at each step should be continuous, 
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and linear forms are preferred to, say, polynomials because they extend 
linearly beyond the domain for which the data are fitted. In addition, the 
selected form should be sufficiently elaborate to meet the accuracy objectives 
of the fitting process; that is, v* should remain within the range of accept- 
able nozzle angles over as much of the acceptable flight envelope of v* as 
is possible; no interior region of this envelope should become unacceptable as 
a result of inaccuracies in the approximating function; and the fitting 
accuracy should be greatest near the boundaries of the acceptable flight 
envelope where the range of acceptable nozzle angles is smallest and control 
power is marginal. 

In the present problem the behavior of v* differs distinctly among the 
regions shown in figure 21(b) but is consistent within any region, so that a 
different function is expected to be required in order to model accurately 

on each region. A suitable form was derived as 

v*(x) = min^l04,max{£Q(x) ,il^(x) ,q(x) ,6°}} (42) 

where q(x) models v* in the region of unconstrained maximum control' power 
and is quadratic in and £(), are both linear in Au and model v* 

in the region of minimum lift margin. Two functions are required in this 
region in order to model the effects of piecewise linear variations in the 
underlying flap schedule. Otherwise, v* is bounded by the minimum and 
maximum nozzle angles. The fitting function is further defined in figure 22(a) 
together with a subdivision of the domain into regions on which v* is given 
by one or another of the functions in equation (42) and which can be compared 
with the subdivision shown in figure 21(b). 

A contour plot of v* is shown in figure 22(b) and can be compared with 
the original function in figure 21(a). 

The approximating function, v*(x), is defined ever 3 rwhere in S both 
inside and outside the acceptable flight envelope ^aF of figure 21, even 
though no data outside that envelope are available since the optimum config- 
uration is undefined for this region in the problem formulation given in the 
earlier section on configuration optimization. It should be noted that trim 
solutions exist in a limited region, outside the acceptable flight 

envelope, even though all such solutions are operationally unacceptable for 
steady flight. For this discussion it is useful to define these flight 
envelopes as volumes in the space of path variables 

•^p(p) = {x; ^ (j), X ^ 

^p(p) = <l>, X E 

The constraint abuse buffer region in which trim solutions exist but are 
unacceptable is then 
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The boundaries of and §)Cy corresponding to v* and F* are 

shown in figure 22(c). It is seen that a buffer region is provided by this 
schedule. This buffer zone is subdivided in figure 22(c) according to the 
constraint parameters which are abused. Along the inner boundary minimum 
speed is limited by lift margin, throttle, and control margin; maximum Au by 
pitch attitude, throttle, and control margin; and minimum by control 

margin and pitch attitude. The outer boundary is reached generally at the 
engine limits. 

In flight, the aircraft may be carried into this buffer region unexpect- 
edly by wind disturbances (gusts, wind shears, and step vertical, or hori- 
zontal winds) which induce abuses of the minimum acceptable speed boundary, 
or more generally, flight envelope abuses. A sudden loss of trim solutions at 
this boundary would create a difficult monitoring problem for the pilot and 
some explicit criterion for defining the configuration schedule in the buffer 
region is desired. The optimization procedure given in the section on configu- 
ration optimization can be extended to the region where no acceptable trim 
solution occurs by relaxing the constraints listed in figure 15, in some 
specified order and amount, until a "least unacceptable" solution occurs 
satisfying the relaxed constraints. The inclusion of such criteria for the 
buffer region in the configuration optimization procedure incorporates the 
design goal that 


ifp(p) = 5^^p(p) 

explicitly in the problem. In the present work, v* was not explicitly 
designed for this region, but an acceptable design occurs due to the ability 
of the piecewise linear fitting function used here to extend in a physically 
reasonable way beyond the domain of the data to be fitted. 

The boundary of has been included in figure 22(c) as a dashed line 

which can be compared with the boundary of achieved by the schedule. The 

two sets differ only along the minimum speed boundary where the schedule loses 

4 knots or less in the descent regime and more in the climb region. 

The locus of operating points in the ^u “ Vp plane for a representative 
AWJSRA STOL approach flight from cruise to glide slope and touchdown is 
included in figure 22(c). Very nearly the entire sequence of descending or 
decelerating segments is within the acceptable envelope ^aF^Ps^ 
designed configuration schedule, but approaches the minimum speed and A^ 
boundaries very closely with a negligible abuse of the minimum speed boundary 
during flare. A significant aspect of STOL approach operations, such as this 
one, is the full exploitation of the low speed and maximum deceleration regions 
of the flight envelope to achieve the steep approach and short landing field 
objectives. 
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COMPLETE TRIM SOLUTIONS 


Solutions for the nominal values of the remaining control variables and 
other configuration-dependent parameters of interest in the control problem 
can now be calculated over the complete flight envelope corresponding to the 
configuration schedule, v*, F*. Contour plots of the nominal values o-f 
throttle, angle of attack, elevator, lift margin, and control power are given 
for the Au - Vg plane at = 1 (fig- 23) and for the Au - plane at 
80 knots airspeed (fig. 24). In flight, values of these controls will be 
adjusted from their nominal values due to random corrective acceleration 
commands, nonstandard parameter values, and inaccuracies in the aircraft 
modeling data used to generate these plots. 

Contour plots in the ^u - Vg plane give the variation of control set- 
tings for all static equilibrium and steady decelerating or accelerating 
operating points. The throttle setting (fig. 23(a)) shows the conventional 
trends in the region of level and ascending/accelerating flight where the 
conventional nozzles-up configuration occurs (cf. fig. 22(b)); that is, 
throttle increases principally with A^ and somewhat with Vg. In the 
region of descending/decelerating flight where nozzle variations occur the 
throttle is more nearly independent of both A^ and Vg out to the edge of the 
envelope. At speeds below 70 knots, where engine power is used to maintain 
lift margin, the throttle is seen to rise rapidly with decreasing speed and 
is nearly independent of A^. Thus, three regions with different throttle 
behavior occur, corresponding to conventional and nonconventional configura- 
tions and to minimum lift margin. 

As in CTOL aircraft, angle of attack principally controls lift coeffi- 
cient and normal acceleration at all flight conditions and can therefore be 
expected to vary directly with A]^, inversely with Vg, and to be independent 
of A^. These trends can be modified by the variation of flap since flap also 
controls lift coefficient, and by the variation of throttle when flaps are 
down since engine power then controls lift coefficient to a significant degree. 
Thus, details of the angle of attack variation seen in figure 23(b) are due in 
part to the selected flap schedule (fig. 17) and to the nozzle schedule at low 
speeds through its effect on the throttle setting. In a manual approach with 
CTOL aircraft the flap history is typically a step function versus speed with 
a corresponding sawtooth angle of attack function. For an automatic approach, 
a smooth flap schedule can be expected to give a nearly constant angle of 
attack. A similar trend occurs with powered-lift aircraft and is extended to 
speeds below that for which maximum flap is reached by raising throttle 
smoothly with speed reductions. Figure 23(b) shows that angle of attack varies 
very little with airspeed or ^u over nearly the entire flight regime; a 
"bowl" at 90 knots and a "peak" at 80 knots are related to corresponding 
variations in the scheduled flap with Vg and A^, respectively, in these two 
regions . 

A consequence of flap usage in CTOL aircraft is that the trim value of 
the regulator control variable a is maintained roughly centered in its 
allowed range; this relationship between scheduling of the configuration 
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control variables and regulator controls also appears in the present maximuin 
control margin design. The nozzle schedule, in that region where the nozzle 
is not at an extreme, maintains the throttle roughly centered in its range, 
and the flap acts similarly on a. 

The trim elevator setting (fig. 23(c)) does not saturate its acceptable 
range anywhere or impose a boundary on the flight envelope. There is little 
variation in trim elevator over the flight envelope. 

Lift margin (fig. 23(d)) increases with speed and is very nearly indepen- 
dent of either or flight-path angle. Along the approach path superposed 

on figure 23(d), this control parameter drops from 2 g at cruise to the 
minimum permitted value for powered lift, 0.4 g, which is reached on the glide 
slope. For the AWJSRA, the lift margin constraint limits the minimum approach 
speed. At higher speeds, excess lift margin can potentially be traded for 
reduced fuel rates by reducing the scheduled flap value. This is based on the 
increase in aerodynamic drag with increased flap setting at higher speeds. 

Fuel reductions are therefore possible in a region with lift margin above some 
comfortable value, say 1 g, (Vg > 100 knots) and scheduled flap settings above 
the minimum (Vg < 157 knots) . 

Control margin for the constant-nozzle regulator is shown in figure 23(c). 
Control margin reaches zero at the minimum and maximum throttle values allowed 
for regulation and is undefined outside that boundary. Good control margin 
(above 0.25 g) is maintained throughout most of the flight envelope. It drops 
below this in the low-speed descending flight regime and in narrow bands around 
the edges of the flight envelope. During the approach path superposed in 
figure 23(e), its value drops from above 0.5 g to 0.12 g at the glide slope 
operating point. Flight is extended into this marginal region with acceptable 
control margin (0.25 g) by using three controls for regulation; otherwise, 
acceptable approach operation with the fixed-nozzle regulator would be limited 
to shallow descents and higher approach speeds enclosed by a 0.25-g contour in 
figure 23(e). Nozzle scheduling to maximize control margin is of course more 
critical in the marginal regime. 

The fixed nozzle mode has low control margin and high control sensitivity 
in the low-speed descending flight regime. This occurs because perturbation 
forces generated by variations in a and 6 ^ are both close to the normal axis 
when nozzles and flaps are down so that there is poor control margin in 
longitudinal acceleration. The fixed throttle mode has more nearly orthogonal 
perturbation forces associated with variations of a and v at high nominal 
nozzle settings and is, therefore, of interest as an alternative central 
regulator mode in this flight regime. To evaluate the potential improvement 
in central regulator mode control margin, data for a throttle schedule which 
maximized the fixed throttle control margin (CMg^^^ 5 ^) was generated by the 
same method previously outlined for the optimum nozzle computation. The 
results showed that the fixed throttle mode had greater control margin over 
the critical STOL regime (roughly, Vg < 100 knots, A^ < 0); the improvement 
was only modest at the glide slope (0.03 g) due to limitations on nozzle 
travel, but exceeded 0.1 g in much of the regime. 
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Contour plots on the “ Ajg plane are given in figure 24 for an 
airspeed of 80 knots. These plots illustrate the variation of nozzle schedule 
and trim controls with A]^ and on various constant speed maneuvers in which 
Ajj differs from 1 g, such as steady turns and flight-path angle changes. The 
flap schedule (eq. (39)) is not included in these plots since it is independent 
of Aj^. 

Figure 24(a) shows some overview characteristics of the configuration 
schedule at this speed. The boundary of trim solutions is irregular and 
defined variously by engine thrust and aerodynamic lift limits. The boundary 
of acceptable flight conditions is defined in different segments by constraints 
on minimum lift margin and minimum or maximum throttle or pitch attitude. The 
intervening buffer zone is hatched to indicate the abused constraint param- 
eters. Operating points for some flight conditions of interest are superposed 
on this figure; the loci of static equilibrium and steady turns for different 
flight-path angles is shown together with a scale of equivalent flight-path 
angles below the abscissa of the figure. A range of angles from -10° to 8° 
is within the acceptable flight envelope at this speed. Flight-path angle is 
changed by transient pullup or pushover maneuvers requiring excursions in 
normal acceleration of 0.1 g to 0.2 g around static equilibrium; maneuvers of 
this type for the AWJSRA approach are shown in the figure and are within the 
acceptable flight envelope. 

The nozzle schedule (fig. 24(b)) varies principally with Au and the 
nozzle is therefore used to change flight-path angle. Lesser variation with 
An occurs so that there is some adjustment of nozzle angle to achieve steady 
turns. There is very little use of the conventional nozzles-up configuration 
at this low speed. 

The corresponding throttle setting (fig. 24(c)) shows little variation 
except near the margins of the flight envelope; hence the throttle control is 
nearly fixed during the maneuvers discussed above and its conventional function 
in controlling A^ is taken over by the nozzles. 

The corresponding trim angle of attack (fig. 24(d)) increases with An 
and, therefore, angle of attack is used to control An during maneuvers 
requiring transient or steady changes in An such as the pullup and pitchover 
maneuvers superposed on figure 24(d), This use of angle of attack is the 
conventional one. The trim elevator setting (fig. 24(e)) decreases with An 
except at positive values of A^ and remains well within the permitted range 
of elevator settings. 

To summarize figure 24, over the maneuver range of interest, the nozzles 
are used to control A^ and angle of attack to control Anj while the 
throttles are nearly constant. This can be compared with conventional con- 
figurations in which nozzles are up and throttle controls A^- More generally, 
the control usage seen here occurs at all lower speeds and at higher speeds in 
the region of the A^^ - plane in which nonconventional nozzle positions 

are scheduled. The region of conventional configuration and throttle usage 
is nearly null at 80 knots but expands at higher speeds. 
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SIMULATION TESTS 


Simulation tests have been carried out in a simplified simulation of the 
system shown in figure 1 (cf. ref, 9). The simulation accurately represents 
the translational degrees of freedom of the aircraft and is suited to testing 
the operation of the trajectory control subsystem in response to low frequency 
inputs such as trajectory commands, Triinmap model errors, and steady winds. 

In this simulation, the combined attitude control subsystem and aircraft 
attitude dynamics are replaced with approximate second-order dynamics for 
attitude response to attitude commands. The estimation subsystem is omitted 
and perfect state estimation is assumed except for a filter to simulate the 
estimation of steady winds. The command generator and path regulator elements 
of the control system are included in the simulation; details of their design 
are outside the scope of this paper but descriptions are given in references 1 
and 2. 


A flow chart of the complete Force Trimmap of the simulated automatic 
control system is shown in figure 25 with relevant equations from the text 
and appendix D noted in each element. The inputs are the reference trajectory 
commands (acr> VC]-) from the command generator, the corrective feedback 
acceleration Aa^- from the path regulator, the wind estimate W;^-, and the 
sideslip command • Additional inputs used are the current flap setting 
(6f) and parameter values (W, x, 6, p) . The vector inputs are assumed 
referenced to runway axes; the acceleration commands are converted to stability 
axis components of the applied specific force (A^^qj^, A^^ , Aj^^) using 

equations (4), (7), (D4) , (D18) , and (D19) , and these are used by the config- 
uration schedule and trim solution algorithm of the basic Trimmap to determine 
^tc 9 ^c^* Finally, the elevator trim setting 6eT is computed from 

the trim solution (appendix A) , and the Euler attitude and angular rate 
commands to the attitude control subsystem Oq, ^^b^ computed from 

various inputs (actp, Yac> ^Vac^ ^Vc> using equations derived in 

appendix D. These inputs are derived in turn from the Trimmap input variables 
as noted in figure 25. 

The approach test path used for the simulation results of this section 
is shown in figure 26 along with the input parameters to the control system's 
command generator which define the trajectory. This trajectory is a simulation 
test path and not one proposed for actual flight. It consists of a sequence 
of legs that correspond to a succession of steady or slowly changing operating 
points. The command generator is required to fill in reasonable maneuvers to 
pass from one steady operating point to the next at the leg junctions. Some 
relatively complex maneuvers involving all three path axes simultaneously are 
included in the test path.' Since nonzero model errors are present in the 
Trimmap due to aerodynamic pr engine model error, programming errors, or 
errors in estimates of airspeed or parameter values, a test is included during 
leg 3 to demonstrate their effects on system performance by off-setting the 
simulated aircraft lift and drag forces by 0.06 g. Finally, steady winds are 
normally present in terminal area operations and wind tests during a level 
turn (legs 4 and 5), a low speed helix (leg 7), and on the glide slope (leg 8) 
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are included to demonstrate the ability of the system to track the four- 
dimensional inertial reference path in the presence of winds and wind shears. 

Simulation results are shown in figures 27(a)-(c), in which the principal 
variables are grouped into acceleration and speed commands (Trimmap inputs) , 
control commands (Trimmap outputs), and trajectory errors (tracking perfor- 
mance) . 


Response to Maneuver Commands 

The initial period up to 55 sec of flight illustrates the maneuvering 
capabilities of the system. Transient acceleration commands to acquire leg 2 
(a 2-axis maneuver to enter a descending turn from leg 1) and leg 3 (a 3-axis 
maneuver to enter a level decelerating straight segment from leg 2) are sup- 
plied by the command generator; the longitudinal acceleration command drops 
without overshoot at the pitchover and then rises without overshoot at the 
pullup to level flight but with some small reversal as the deceleration of 
leg 3 is initiated. The normal acceleration command drops transiently for the 
pitchover and then rises without overshoot to its steady-state value for the 
turn; the converse transient occurs at the exit from leg 2. 

The applied specific force time histories in figure 27(a) include traces 
for the total command, which includes the feedback term, and the actual air- 
craft applied specific force. Since no external disturbances are present 
during this period, transient differences from the nominal specific force 
commands measure the ability of the system to follow maneuver commands. It 
is desired that these response errors be minimized in order that the regulator 
work chiefly on the effects of external disturbances to the system. The 
response errors depend on both the commanded maneuver defined by the command 
generator and on the acceleration tracking accuracy of the Trimmap-aircraf t 
combination. In general, nonzero maneuver command rates (principally jerk) 
require nonzero control rates and these result in trajectory errors and feed- 
back activity due to control lags; that is, control rates and tracking errors 
increase with the commanded maneuver activity. Thus, design goals for the 
command generator include minimizing the maneuvering required to transition 
between any steady flight conditions specified by the coarse input trajectory 
commands to the system, and to limit maneuvering commands in accordance with 
the acceleration tracking capability of the Trimmap-aircraf t combination. 

The activity of the scheduled controls for the first 55 sec shows only 
minor nozzle use which is closely related to the change in normal-specific 
force for the turn, while the regulator controls show throttle variations 
which follow the change in longitudinal specific force with flight-path angle. 
Only minor angle of attack variations are used so that pitch attitude simply 
follows flight-path angle. The roll command is seen to vary smoothly without 
overshoot at turn entry and exit; this behavior reflects the success of the 
command generator design for the lateral axis. 

These maneuvers occur in the CTOL flight regime and begin and end with 
the AWJSRA in its conventional configuration; the control activity seen here 
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differs from conventional control in that nozzle activity replaces angle of 
attack for the turn execution. If operationally desirable, CTOL-like maneuver 
control activity can be imposed on the Trimmap by requiring the use of minimum 
nozzle setting as much as possible in the configuration schedule. 

The control activity used by the Trimmap to execute the commanded transi- 
tion maneuvers in figure 27(b) shows desired characteristics; that is, no 
unnecessary overshoot or reversal. More generally, control activity, as mea- 
sured by control rates required for maneuvering flight, is given by 


u 

c 



A 

3Au ^ 


3u* r 

3A„ - 


[VxU*(x,p) ]x^ 


(45) 


where 


— T 

u — (6^,v,ot,6^) 

This expression contains the gradients of the configuration schedule and 
associated trim solutions of the Trimmap, u*, and the maneuver command rates 
of the command generator, x. Thus, control activity and its minimization 
for the automatic control system depends on both of these elements. 


The control activity seen in figure 27(b) involves simultaneous use of 
roll, angle of attack (or elevator), throttle, and nozzle controls during 
maneuvers requiring accelerations along the three path axes simultaneously; 
that is, the Trimmap automatically coordinates all of the AWJSRA’s controls 
with the commanded maneuver for any maneuver in the output set of the command 
generator. In contrast, manual flight operations are designed to minimize 
pilot workload, and a significant problem for STOL aircraft is to limit work- 
load requirements within the capabilities inherent in manual flight control; 
this generally implies control usage limited to two controls at a time with 
sequential step changes of the redundant configuration controls, and maneuver- 
ing carried out sequentially along one axis at a time (refs. 6, 7, and 8). 

For the present automatic system, these operational restrictions were not 
imposed because the workload capacity of the digital flight-control system, 
as measured by the number of input and output channels, data memory capacity, 
and real-time processing speed, exceed those inherent in direct manual flight 
control, but the objective remains to minimize control activity. Restrictions 
on the operational use of this general capability of the automatic control 
system fbr ease of monitoring requirements and pilot takeover were not con- 
sidered in this study. Restrictions on maneuver complexity are readily 
imposed as restrictions on the coarse input trajectory commands. 


Effects of Steady Model Errors 

The actual acceleration can differ from the commanded acceleration both 
transiently and in steady state due to errors in the model of aircraft forces 
used by the Trimmap to compute control settings. These errors can arise from 
several sources (aerodynamic and engine force model errors, state and parameter 
estimation errors, programming errors, neglected dynamic effects) and are 
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normally present during flight. To illustrate the effects of steady model 
errors, the normal and longitudinal applied specific forces of the aircraft 
were offset by 0.06 g during the time interval [55, 80]. These errors are 
compensated by feedback as follows: in steady state, position standoff errors 

occur (fig. 27(c)) of sufficient size to generate a steady feedback accelera- 
tion command, Aa, which exactly compensates the model error. This effect is 
seen in figure 27(a), where the total specific force commands (Au^j > rise 

to steady values different from the nominal specific force commands, but the 
actual specific force (A^, An) is held at the desired nominal value. The 
required steady corrective acceleration command is generated by steady offsets 
of the regulator controls, (a, 6^) from their expected value as given by the 
configuration schedule and available model data. This is seen in figure 27(b) 
where throttle is driven nearly to its acceptable limit to compensate for the 
longitudinal error which, in size, is a large percentage of the engine’s 
acceleration capacity. The angle of attack compensates for the normal accel- 
eration error but in this case only a very small offset occurs because the 
model error is a small percentage of the wing’s lift force capacity. 

In general, if is the Trimmap model error expressed as an accelera- 

tion error vector, 

= (Au,An)^ 

then a steady position standoff error proportional to the acceleration error 
develops asymptotically 

AR-^[K-l]£a (^6) 

where [K;^] is the position error feedback gain in the path error regulator. 
This induces a steady corrective acceleration that exactly compensates the 
model error 


Aa-^-Eg (47) 

but requires a steady offset of the regulator controls from their values 
predicted by the design data 


Au^ -H^-[VaU^(x,p)]e^ (48) 

Here [VaU*] is the Jacobian of the regulator controls with respect to accelera- 
tion. The derivation of these results is omitted and rests partly on details 
of the path error regulator element in figure 1. The basic dependence of 
path tracking performance on the modeling accuracy of the trim solution algo- 
rithm and the related trade-off between performance and model accuracy are 
evident in equation (46) . Steady feedback signals measure the modeling error 
(eq. (47)), and an analysis of this signal from flight records would indicate 
the size and nature of such errors. In addition, these errors would saturate 
the authority limits of the feedback if sufficiently large. The system 
compensates for model errors by adjusting the regulator controls (eq. (48)) 
and may thereby result in reduced control and lift margins. 
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The path tracking errors and steady feedback offsets noted in equa- 
tions (46) and (47) can be removed or reduced by adding integral feedback 
of the acceleration error to the Trimmap of figure 25; this can be implemented 
as an additional corrective acceleration, 6^, in the total acceleration command 

actp = acp + Aap + Sap (49a) 

where, for each path axis, the correction has the form 


6a 


max*|-6a^,min|K 


(ac - a)dt 



(49b) 


Here, 6a^ is the authority limit of the integral feedback, a is the estimated 
acceleration, and K the feedback gain. This model error compensation device 
adjusts controls in accordance with equation (48) , but accommodates model 
errors |e^| up to the amount without steady-state position error or 

feedback offset. Equations (46) and (47) then define the effects of steady 
Trimmap inaccuracies in excess of A compensator of the type (eqs. 49)) 

is not included in the present simulation study but will be included in a 
further study of the system. 


Accuracy requirements for the Trimmap model are of interest because of 
the relation between model accuracy and Trimmap computational requirements. 

The Trimmap model is derived from an estimate of the aircraft and contains 
errors due to approximations in representing the estimated aircraft and the 
uncertainties of the estimated model. With the use of integral feedback, 

(eqs. 49)), approximation errors up to a fraction of 6a^ can be permitted 
while the remaining compensator authority is applied to the basic uncertainties 
in estimating the actual aircraft. 


Response to Winds 

A steady x-axis wind of 15 knots is imposed during a 180° turn in the 
time interval [100, 155] and, with respect to the aircraft, acts variously 
as a tailwind, crosswind, and headwind during the turn. The control system 
is expected to track a circular path, as seen from the ground, and maintain 
the desired ground speed. 

The appearance of steady wind in the Trimmap calculations is seen in 
figure 25; the velocity vector with respect to the air mass is used to define 
the path axes in which the components of the nominal and total acceleration 
commands are given, and to define the equivalent airspeed used in the config- 
uration schedule and trim solution algorithm. In the simulation results in 
figure 27(a), rises and falls during the turn due to a geometric effect; 

the required air velocity vector is not perpendicular to the radius_ 

vector of the inertial circle to be tracked and some acceleration along 
is required to supply the centripetal turn acceleration.^ These effects of 
winds on the inputs to the basic Trimmap can be derived; the velocities of 
interest are related by 
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= Vi - w 


(50) 


The effect of the estimated wind on the spherical coordinates of is then 


AV. = -u • W 
A — 


Ai|^^ = -u ® . k 


Ay = YtH 


W 

Vt 


A 


(51) 


J 


Here, relatively small winds are assumed 

lw| 


« 1 


and the effect of steady vertical wind, W^, on :^light-path angle iS' neglected 
since Wz normally negligible. _During the turn, is time^varying as 

Vj changes direction. Otherwise, is constant when Vj and W are constant 

vectors. The corresponding variations in Trinmiap inputs are derived from 
equations (6) and (51) as 


Ax = (AVg, AA^, AA]y^) 


AVe = AVa /p/Po 


AA^ = a • Au = a • mAij^y + Ay 


AAj^ = cos (j)^a • An + sin cf)^a • Am = 0 




(52) 


As seen in equations (52) the centripetal turn acceleration a • m, which is 
of the order of 0.4 g typically, will affect the commanded acceleration along 
the aircraft stability axis u^ proportionally to AijjyCt) . The effect on 
normal acceleration command is generally negligible compared to 1 g, by 
several orders of magnitude. 

Because the inputs to the configuration schedule are adjusted for 

winds, the configuration controls will generally be adjusted automatically 
to maintain the designed maximum control margin configuration in the presence 
of steady winds. This is illustrated in figure 27(b), where the flap command 
varies in relation to the equivalent airspeed command although there is* 

no adjustment of nozzle from its minimum value at this example flight condi- 
tion. The regulator controls also vary in response to the effects of winds 
on x^, as seen in the throttle variation which follows the time history of 
Au^. In the present example, note that there is little effect on angle of 
attack since the effect of speed variations on lift is taken up principally 
by the scheduled flap adjustments. The general effects of steady winds on 
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adjustments of the controls from their anticipated no-wind values can be 
derived from the configuration schedule and associated trim solutions 
(figs. 17, 22(b), and 23) and can be represented as 

Au = [V^u*(x,p) ]Ax (53) 

where Ax is given by equations (52) and [V^^u*] is the Jacobian of the basic 
Trimmap function or algorithm. At flight conditions other than the test case 
of figure 27, the nozzle may be adjusted in place of throttle (when v* is 
not at a nozzle limit) and, similarly, a may be adjusted in place of flap 
(when F* is at a flap, limit); in all cases, however, the Trimmap maintains 
the optimum configuration automatically. 

The trajectory tracking performance illustrated in figure 27(c) shows 
good tracking of the commanded velocity vector of the trajectory (Vj^, , Yc) 

during the wind test with only minor transient errors due to transient maneu- 
vering at turn entry and exit and due to wind estimation lags at wind onset 
and offset. 


Pitchover and Configuration Transition 

During the period, t ^ [160, 194], the aircraft exits a decelerating turn 
(leg 5) onto a steep descent (leg 6) . The longitudinal acceleration command 
(fig, 27(a)) reaches its descent value without overshoot while the normal 
acceleration transient for the pitchover uses the maximum excursion, 0.125 g, 
permitted for flight-path angle changes in the command generator. 

Large changes in the configuration controls occur during the pitchover 
(fig. 27(b)); flap rises to 55° and nozzle to over 80°. The scheduled flap 
variation for this maneuver is due to the speed reduction and is independent 
of the flight-path angle change (fig. 17) while the nozzle change is due 
principally to the flight-path angle change (fig. 22(b)). The regulator 
controls, a, 6|_, are used transiently for the pitchover and are adjusted 
somewhat for the descent but remain roughly in the center of their ranges. 

Good maneuvering commands and execution is evidenced in the roll angle and 
in the velocity coordinates and path errors of figure 27(c)). 

The aircraft begins the STOL phase of its approach at the pitchover, 
characterized by low speeds and steep descents near the minimum boundary 

of the acceptable flight envelope and by the powered-lift, nozzle-down config- 
uration. Control characteristics such as control margin (fig- 23(e)) and 
sensitivity also change rapidly for the worse with the configuration and 
speed changes and will remain marginal until landing. 

This interval of the test path can be compared with the corresponding 
transition maneuver for manual approach described in reference 6, in which 
flap is first reduced stepwise with speed in level flight and then a step 
change in nozzle to its approach setting occurs simultaneously with pitchover 
to the descent path. In the present automatic approach, a short transition 
period again occurs, separating longer conventional and STOL phases of the 
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approach, and the flap and nozzle variations show the same relationship to 
speed and equivalent flight-path angle. However, the automatic transition 
shown here does not attempt to minimize pilot workload, so the sequential and 
descrete steps of the manual transition are now carried out simultaneously 
and continuously. The automatic system attempts to minimize control activity; 
that is, overshoot of steady-state value or unnecessary reversal is avoided in 
the path and control variables (A^, V, v, a, 6^) . 


Descending Helical Turn with Winds 

Leg 7 (t ^ [194, 294]) is a low-speed, steeply descending and decelerating 
turn in the presence of a steady 10- knot y-axis wind whose onset occurs at 
the quarter turn (t = 205) . The maneuver execution is similar to earlier 
results except for wind effects. The 10-knot wind is a crosswind relative 
to the runway center plane and to approach paths confined to the centerplane, 
and is not difficult to control as a crosswind relative to the aircraft. 
However, during the 360° helical turn it acts variously as headwind, crosswind, 
and tailwind and produces an excursion of the equivalent airspeed Vg down 
to 58 knots in that quadrant of the helix in which it is a tailwind 
(fig. 27(a)). This excursion abuses the minimum speed boundary by 6 knots and 
results in abuse of the lift margin and maximum throttle constraints. In 
addition, the longitudinal acceleration is affected by the wind 

through both the turn acceleration and flight-path angle terms in equa- 
tions (51) and (52) with excursions of Au down to -0.18 g. This excursion 
abuses the minimum A^ boundary of the acceptable flight envelope during 
about half the turn and results in very low control margin. These envelope 
abuses are also shown in figure 28. The minimum approach speed, » 

is a familiar flight safety parameter for CTOL aircraft; in addition, STOL 
aircraft can approach at steep glide angles near the minimum longitudinal 
applied specific force boundary so that Au^^^^(p) is also a significant safety 
parameter. 

The control activity (fig. 27(b)) shows that the effects of wind on lift 
are balanced at first by flap adjustments until it saturates at its maximum; 
then by a increases, and finally, at the lowest speeds, large throttle 
increases are used to supply powered lift and a modest abuse of the maximum 
allowed steady throttle (25.6°) occurs. The nozzle variation is related to 
the effect of wind on rises to values above 90° where there is 

little speed control margin remaining in the regulator. The ability of the 
Trimmap and autopilot to survive envelope abuse events and to continue solving 
the trim equations outside the boundary of the acceptable flight envelope has 
been demonstrated here. This capability is an important flight-safety aspect 
of the design and, although the configuration schedule design procedure 
formulated in the section on configuration optimization does not supply design 
criteria for the envelope abuse regime, such criteria can be incorporated in 
a general design methodology for redundantly controlled aircraft. 

In the present simulation tests, the system is carried outside its 
acceptable envelope while attempting to track the commanded four-dimensional 
inertial path in the presence of winds. In an operational environment, a 
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four “dimensional guidance system would be additional element of the automatic 
control system not present in the system of figure 1 which would synthesize 
trajectory commands using wind information so that the aircraft is not forced 
to abuse its acceptable operational envelope in order to track the inertial 
trajectory (e.g., refs. 11 and 12). A helical path in the wind environment 
of the present simulation test can be flown at a shallower descent angle and 
higher speed without abusing the flight envelope, although the susceptibility 
of four-dimensional helical turns to wind effects on control activity and con- 
trol characteristics makes it a doubtful operational maneuver. 


Glide Slope with Wind Shear 

The glide slope, (t ^ [294, 381]), is a steep descent at minimum speed, 
which is disturbed in this test by both a crosswind (10 knots) and a headwind 
(15 knots). The headwind is initiated at t = 315 and then both winds 
decline linearly to zero at t = 360, prior to flare. This shear is modest 
relative to values encountered in practice but suffices to illustrate the 
effects of steady shear. 

The crosswind requires only a change of aircraft heading ip from the 
desired inertial direction this correction is supplied automatically 

by the Trimmap, as seen by the small difference between ipy and iIj in fig-- 
ure 27 (c) . 

The headwind produces variations of the nominal commands; both and 

increase temporarily in accordance with equation (52) and the wind 
time history and also transient feedback components occur in Au^ and Aj^^ at 
the onset of the headwind (fig. 27(a)). The feedback components arise both 
from lags in the control response to the wind variation, and from errors in 
estimating the wind. Wind estimation errors enter the Trimmap calculations 
of lift and drag forces as airspeed errors and, hence, produce the same effect 
as other modeling errors (eqs. (46)-(49)); that is, if <5V is the steady 
airspeed estimation error then the corresponding steady acceleration error is 



In the glide slope test, steady wind estimation errors result from wind shear, 
and the consequent steady position offsets below the reference glidepath, 

^Rn’ seen in figure 27(c). Larger steady shear values would produce 

proportionately larger estimation errors, 6V, and corresponding position 
offsets . 

The control activity on the glide slope (fig. 27(b)) shows some nozzle 
adjustment in response to the wind effect on A^. In manual operations 
(ref. 6) the nozzle is adjusted on the glide slope for winds; here, the 
automatic configuration schedule carries out a similar adjustment and can 
vary the adjustments with wind changes due to shear. The configuration 
schedule also attempts to reduce flap as airspeed increases, but is prevented 
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from doing so by a flap retraction stop installed in the aircraft (and simu- 
lated here) that prevents flap retraction by the automatic control system 
below a detent, once that detent is passed during an approach. The detents 
occur at flap settings 30°, 50°, and 65°. The retraction-stop event has the 
same effect on the trim solution calculation as a flap rate saturation event; 
that is, 6f^ and 6f are unequal and can remain unequal persistently in the 
case of steady winds. The calculation of the regulator controls, (a, 6t) > 
should therefore be based on the actual flap in order to generate the com- 
manded acceleration accurately, as is done in the present Trimmap design. 

Both regulator controls show variations on the glide slope (fig. 27(b)) 
due to winds; a complex throttle transient occurs at wind onset which in part 
executes the feedback transient in at wind onset and, in part, the 

nominal variation (see fig. 27(a)). In the STOL configuration, engine 

power is used to generate both normal and longitudinal acceleration and, 
hence, even if each input signal is a well-behaved transient, their combined 
effect on the throttle command can be twice as complex and require a higher 
frequency, as seen here. The angle of attack variation follows the airspeed 
change and also is used by the regulator function of the Trimmap to supply 
the corrective feedback portion of the combined effect is a rapid drop 

in cx and pitch attitude by 5° at wind onset. 

Several overview figures of the Trimmap inputs and outputs on the test 
path are shown in figures 28(a) and (b) . The trajectory variables 
{V£(t), Au(t), t ^ [0, 384]} are seen in figure 28(a), and a dashed trace 
of the trajectory in the absence of winds is included for comparison. It is 
apparent that the winds during the level turn, helical descent, and glide 
slope legs can have a pronounced effect on these Trimmap and configuration 
schedule inputs while tracking a four-dimensional approach path. The low- 
speed steep helical descent is especially susceptible to envelope abuse due 
to winds both in minimum speed and minimum equivalent flight-path angle. 

The locus of configuration commands (fig. 28(b)) during the approach is 
relatively simple and migrates from the lower left corner of the figure 
(CTOL configuration) to the upper right (STOL configuration). Of the remain- 
ing corners, the upper left does not occur in the configuration schedule and 
the lower right corresponds to nominal flight conditions which are little used 
in passenger operations. However, these regions would come into greater use 
if the nozzle were used for regulation. 

It should be noted that, although some feedback activity and exercise of 
the regulation function occurs in the simulation results presented in this 
paper, these tests primarily demonstrate the functioning of the automatic 
configuration schedule and autopilot system in response to low-frequency inputs 
that excite only a minimum of feedback activity. Realistic examination of the 
regulator behavior and response to higher frequency maneuvers and feedback 
commands and to higher frequency disturbances and modeling errors requires 
a more accurate six degrees-of-f reedom simulation of the aircraft attitude and 
engine dynamics than was used here; such an examination will be the subject of 
future simulation tests . 
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CONCLUSIONS 


The design criteria for a Trimmap for the application of an advanced auto- 
matic control system structure (ref. 1 and fig. 1) to the augmentor wing jet 
STOL research aircraft (AWJSRA) has been described. This structure accommodates 
the coupling of redundantly controlled aircraft to an ATC system having a 
general set of approach trajectories at its command. The Trimmap is the 
central element of an acceleration controller structure capable of executing 
trajectory commands. Acceleration commands are constructed as the sum of the 
acceleration on the commanded trajectory and a feedback term to correct path 
errors; the Trimmap then calculates the control settings required to generate 
these accelerations at the current state and parameter values, using a model 
of the aircraft forces. The principal problems to be resolved in the Trimmap 
for this and other STOL aircraft are control redundancy and the nonlinearity 
of the trim equations. 

A structure for resolving redundancy was proposed (fig. 10). It contains 
two principal elements: first, a configuration schedule drives the redundant 

controls (flap and nozzle) as a function of the nominal trajectory command and 
is configured off-line. The second element is a trim solution algorithm that 
carries out the regulation function by solving the two trim equations for the 
remaining controls (angle of attack and throttle) on the basis of the total 
acceleration command. It was found that the control margin available for reg- 
ulation is maximized by using as many of the controls for this purpose as have 
sufficient bandwidth. The trim solution algorithm was therefore structured to 
default to nozzle use when either of the conventional controls saturate (fig. 14). 

The configuration can be selected optimally at any flight condition 
within a range limited by the operational constraints. For the AWJSRA, the 
flap and nozzle were selected to maximize lift margin and regulator control 
margin. This was done because these margins are strongly dependent on the 
configuration controls and in order that the automatic system retain the 
aircraft’s full capabilities for low-speed, steeply descending flight. It can 
be concluded that configuration management is a key design issue for redun- 
dantly controlled aircraft and governs the operational envelope of the auto- 
matic control system as well as the behavior of configuration-dependent 
parameters over that envelope. 

Aircraft nonlinearity is treated in a simple way by the computational 
methods described here for solution of the trim equations (appendix B) and 
for optimization of the configuration schedule (appendix C) . Tabulated data 
is used so that the relevant nonlinearities are approximated as continuous 
piecewise linear functions with accuracy governed by the number of pieces 
used. In addition, the design domains for all independent variables in the 
problem are bounded and known in each design problem. Thus, solutions of 
the nonlinear trim equations can be computed routinely and reliably by 
representing the model as piecewise linear, using a finite number of pieces, 
and then searching the pieces exhaustively for solutions. This approach, 
however, requires specially devised algorithms for each application and 
more general algorithms based on Newtonian iteration are attractive 
alternatives but of unknown convergence. 
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Similarly, th^ configuration optimization procedure is simply one of enumerat- 
ing trim solutions and constraint parameters on a set of points which span the 
design domains of all variables defining configuration and flight condition, 
and selecting the optimum configuration at each flight condition from the 
solutions obtained* 

The automatic control system, including the Tr immap described in this 
report, was subjected to simulation tests using a four-dimensional (4-D) STOL 
approach path (fig. 26). 

These investigations demonstrate the system response to low-frequency 
inputs, including maneuver commands, Trimmap model errors, and steady winds, 
but system behavior in response to higher frequency disturbances, such as 
gusts and transient model errors, are not covered in these tests. Simulation 
results are shown in figure 27 from which a number of conclusions were drawn. 

1. The proposed Trimmap structure provides automatic coordination of 
all aircraft controls during transition maneuvering and does so for any 
maneuver in a general output set permitted to the command generator and the 
underlying ATC system. Control rates and system performance during transi- 
tions depend on both the Trimmap configuration schedule and the command 
generator. The automatic system has a higher workload capacity than is 
acceptable in manually controlled flight, and the present concepts exploit 
this higher capacity to achieve the generalized maneuver capability. 

2. Steady Trimmap errors in modeling the aircraft forces are continu- 
ously present in flight and can arise from aircraft model data inaccuracies, 
state or parameter estimation errors, and software inaccuracies. The system 
feedback adjusts the regulator controls from their nominal values in order to 
compensate for such errors, and results in a steady position error in 
proportion to the accuracy of the model used by the trim solution algorithm 
(eqs. (46)-(49)). It was concluded that integral feedback of the acceleration 
error (eq. 49)) should be added to the Trimmap in order to relieve the 
regulator of using its authority for steady model error compensation and to 
reduce the sensitivity of system tracking performance to model accuracy. 

3. Steady winds are normally present throughout the approach and can 
significantly affect the Trimmap inputs compared to their values in the 
absence of winds (eqs. (51) and (52) and fig. 28(a)). The Trimmap adjusts 
all controls for the wind effects to maintain the maximum control margin 
configuration automatically while tracking the four-dimensional inertial 
path. Steady wind estimation errors arise during wind shear conditions and 
produce the model error effects described above. 

4. The STOL approach path operates near the boundaries of the acceptable 

flight regime both in minimum speed and in minimum longitudinal applied 
specific force in order to achieve the short field length capability. The 
minimum flight-path angle boundary corresponding to therefore. 


46 



a significant safety parameter for STOL aircraft in addition to the minimum 
speed, Abuses of both boundaries can arise due to winds, and the 
capability of the Trimmap to operate during such abuses was demonstrated. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif, 94035, April 12, 1978 
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APPENDIX A 


AWJSRA ENGINE AND AERODYNAMIC FORCE MODELS AND ELEVATOR TRIM FIRTCTION 


The engine and aerodynamic force models used in the calculations of this 
paper are described in detail. Both graphs and tabulated data are given for 
the various functions required for the models (Tjj, Tq, mg, Cp, 
addition, elevator trim values and engine fuel flow data used in the Force 
Trimmap are ir^luded. 

Function evaluation is carried out by interpolation of the tabulated data 
using simple, invertible interpolation formulas given in appendix B. To 
minimize storage and computation time requirements, table sizes have been 
approximately minimized within the accuracy requirement that the original 
data, based on theory, wind-tunnel tests, and flight tests, be reproduced by 
the force model to an accuracy of 0.01 g. Other classes of data-fitting 
functions (spline functions and orthogonal polynomials) were investigated 
for computational efficiency; however, they were difficult to invert and 
required as many parameters as the tabulated data method in order to obtain 
^he required accuracy. In addition, tabulated models permit straightforward 
control over the maximum error in the model while analytical function models, 
whose parameters are computed via least squares techniques, are designed to 
minimize average errors. The table interpolation method is applied consis- 
tently throughout this paper and was found to be simple and routinely appli- 
cable in all phases of the work. 


Engine Model 

The AWJSRA is powered by two twin-spool jet engines (Rolls Royce Spey 
MK 801-SF engines). Some air is ducted off the first compressor and exhausted 
through the augmentor flaps to obtain lift augmentation; the air is also 
ducted to the ailerons for boundary-layer control. The remaining air passes 
through the engine normally but is exhausted through movable nozzles to 
provide a jet thrust vector which can be varied in direction over 90®. The 
engine nozzles and forces are illustrated in figure 29. 


The direct engine forces are the thrust from the hot exhaust and the 
ram drag due to momentum rate of the inlet air. These forces, referred to 
stability axes, are 


( cos(v + a)\ /“^\ 

)+Va( ) 

sin(v + a / ^ 0 ' 


(Al) 


In addition, the effect of the ducted air used for lift augmentation 
is included in the aerodynamic description through the dependence of lift and 
drag coefficients on the "cold thrust coefficient" 
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(A2) 


C 



where Tq is the equivalent thrust obtained from isentropic expansion of the 
compressed ducted air to ambient atmospheric conditions. 


Three parameters (Tjj, T^,, m£) in equations (Al) and (A2) remain to be 
defined. These are functions of airspeed, atmospheric pressure and tempera- 
ture, and engine rpm. Plots and tables are given in figure 30 and in table 1 
for the "corrected” engine parameters 


Tr(N^//t,V) T^(Nr//t,V) /t mgCNR/Zr) 

~6 ’ ' ”6 ’ 6 


where 


t,6 


temperature and pressure ratios, — 

^SSL 


P 

P 


1 

SSL 


PssL^'^SSL = standard mid-latitude sea level atmospheric pressure, temperature 
(10.13 N/cm2, 288.15° K) 


Pl,Ti = engine inlet stagnation conditions 

These data derived from engine test stand data from the manufacturer and from 
company reports from the aircraft modification contractor, include the effects 
of power bleed and ducting losses. The ratios t, 6 are convenient similar- 
ity parameters introduced so that engine thrust and mass flow at all atmo- 
spheric conditions can be calculated from those at the reference conditions, 

T = 6=1. The independent variable, N^, is the cockpit tachometer dial read- 
ing ("power setting") which, for this installation, is rpm/121.3. Engine 
operational limits on power setting are 


Idle power 

II 

H 

59.3 




Maximum continuous power 

%IC " 

95.8 


> 

(A3) 

Normal T.O. power 

Nto 

98.5 




Emergency power 

Ne = 

101.4 

J 




Such limits occur for most jet engine installations; they are selected to 
yield specified time intervals between required major engine overhauls and 
are marked by detents on the throttle controls. 

The throttle lever angle, which is used to control engine output, is 
directly related to power setting for the AWJSRA (fig. 30). Engine fuel flow 
is also given in figure 30 and table 1. 
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The domain over which these data define the engine model 

0 V < 250 knots 
0 < N„//r < 103.5 

suffices to cover the range of airspeeds of interest for the AWJSRA and opera- 
tion at or below maximum continuous power for temperatures above -12° C. 

The maximum steady thrust available in flight occurs at maximum continuous 
power, it decreases with altitude and temperature and increases with 

airspeed, as shown in figure 31. The maximum cold thrust coefficient also 
occurs at It varies principally with airspeed (fig. 31(c)) and is 

independent of pressure altitude. As seen, Cj is restricted to low values at 
cruise speeds (160 knots) and there is little effect of engine power on aero- 
dynamic coefficients at higher speeds. Much higher values of Cj are avail- 
able at landing speeds (65 knots) and at low speeds generally as desired for 
the lift augmentation system design. 


Aerodynamic Forces 

The aerodynamic forces in stability axes are written 



(A4) 


where C^, Cp are the usual lift and drag coefficients and is the wing 
area, 80.4 m^. 

The aerodynamic description of the AWJSRA used here is derived from 
reference 4, where wing-body lift and drag coefficients are tabulated as 
functions of three variables, {a, 6£, Cj}, over the domain 


a e [-10. 5°, 27.5°] ^ 

6£ e [5.6°,72°] 


(A5) 


The data of reference 4 were based largely on wind-tunnel data, but more 
accurate data have since been obtained from flight tests (ref. 6) and are 
included in the model of this paper. In addition, the coefficients have been 
corrected for the symmetric deployment of aileron with flap that is built into 
the flap mechanism, and for the tail lift which is required to trim the 
aircraft. Ground effects, which occur at altitudes below 16 m are not 
included in the model. 
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(A6) 


The tail lift correction is computed from 

^ + ^lClwb - 

Here, wing-body pitching moment, c is the mean aerodynamic 

cord, and a.p, a^, and aj) are the moment arms for tail lift and wing -body 
lift and drag forces (ref. 4). The tail lift correction contains only func- 
tions of {a, 6£, Cj} and can be included in the lift coefficient without 
increasing the number of independent variables. Variations in c.g. location 
with aircraft weight and the pitching moment due to engine thrust have only 
negligible effect on the tail lift and are neglected. 

The resulting aerodynamic model is tabulated in table 2 and pictured in 
figure 32 as lift-drag polars for four flap settings, with lines of constant 
a and Cj mapped on each plot. During approach and landing the flap setting 
progresses slowly from 5.6° at cruise to 65° at landing. At cruise, engine 
power is seen to have little effect on lift and reduces the aerodynamic drag 
by exhausting pressurized air rearward through the flap. In any case, there 
is little available range of Cj at cruise speeds and correspondingly little 
effect of engine power so that the resulting model at cruise closely resembles 
that of a conventional aircraft and has a maximum of about 2. At low 

speeds and maximum flap the available range of Cj is much greater and its 
effects on lift coefficient is pronounced with maximum in the range of 4 

to 5 . 


The purpose of a conventional flap is to increase the maximum Cl avail- 
able from the wing in order to permit lower landing speeds. The special flap 
of the augmentor wing has the same purpose but obtains a higher Cl^^^^x than a 
conventional flap in order to obtain sufficiently low landing speeds for short 
field length landings. Thus, the polar plots in figure 32 although roughly 
similar to those for conventional aircraft, require an additional parameter, 
Cj, to account for engine power effects; they also exhibit large nonlinear 
variations with flap setting. As a result, the data required to define the 
aerodynamic force model is greatly increased compared to the single lift-drag 
polar required for CTOL aircraft. 

Elevator Trimmap 

In the analysis of the trim equation in this paper the control variables 
are taken as 


u = (6f,6t,a,v) 

of which (6£,6 |-,v) are directly controlled by the autopilot through servos, 
and a is controlled by commands to the elevator servo. Therefore the rela- 
tionship 6e(u) is necessary for the implementation of the autopilot and con- 
straints on elevator usage will correspondingly limit the usable range of u. 

The elevator affects the aircraft forces and moments through the tail 
lift, which is modeled by the linear relation 
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(A7) 


L'p = 


+ L. 6e 
6e 


with elevator travel limited to 


6e e [-25%15°] 


and 


(A8) 


6e 


- QS C. 
^ w L 


6e 


Elevator stall can be neglected since it has not been observed for angles of 
attack below wing stall. In addition, reduction of dynamic pressure at the 
tail due to the wing wake presence is also negligible because of the T-tail 
design of the aircraft. L-p is the tail lift at zero elevator angle and can 
be written as 


LTq QStCLx 

a 

where 


Aotf — ot 


c + L*j> 


T 


tail downwash angle 
tail incidence 

angle of attack for which tail lift is zero 

9Cl^ 

8a 

The downwash angle is generated from a complex model (ref. 4) whose 
independent variables are a-p, Cj}. Here, a-j- is the wing angle of attack 

delayed by the time it takes for air to reach the tail from the wing: 

= a(t - At) 



At 



where a^ is the distance between wing and tail center of pressure. For 
steady flight a-p and a are equal. Values of the remaining parameters of the 
tail lift model, • • -j otop, are given in reference 4. 
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The elevator is used to control the pitching moment balance, given by 


■*■ ^t(LTo ■*■ (A9) 

Here, C.^ are dimensionless engine forces 


C 


DR 


®E^a 

QSw 



and a-j-, ag, a^ are the moment arms of the engine and tail lift forces. 

These moment arms are computed from a nominal center of gravity location. 
Changes in this location with aircraft weight are small and can be neglected 
in computing the tail lift and elevator setting for trim. Equation (A9) can 
be solved for 6e and the resulting terms arranged in three parts to separate 
the effects of dynamic reactions, engine forces, and aerodynamic forces and 
moments . 


6e = 

6ej) = 
6eE = 

6e. = 


Seo + <5eg + 6e^ 

(qly + prdx - Iz) + (p2 - r 2 )I^ 2 > 

Sw 

"V^%B ■ ^D*^DyB) CL^gS^a^ 




y (AlO) 


Cl 


Aot-a 


6e 




In the present analysis the elevator setting for steady or very nearly 
steady flight conditions along the nominal path is of interest so that the 
dynamic term can be neglected and then the elevator setting for trim is 


6e^ = 6e^(a > >Cj) + 6 eg (All) 

The principal term in equation (All) is which is a function of only 

{a,6f,Cj} and can be calculated and tabulated over the same grid for which 
the lift and drag coefficients are given. The results are shown in figure 33 
and table 3. Ground effects (which are of a significant size) are omitted 
from 5eA* 


The physical limits of elevator travel (eq. (A8)) limit the range of II 
for which trim is possible. Under dynamic conditions, the elevator is also 
used transiently to control attitude and therefore the allowed trim settings 
must be further restricted within equation (A8) to provide a margin for the 
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attitude control function. For the AWJSRA design this margin was taken simply 
as 8° and it was found that the constraints on elevator usage had almost no 
influence in establishing the boundaries of the flight envelope or the range 
of acceptable configurations at any flight condition within the envelope. 
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APPENDIX B 


SOLUTION OF THE TRIM EQUATIONS 


Solution of the aircraft trim equations for two control variables is 
required extensively in computing the configuration schedule and is proposed 
as part of the Trimmap element of the flight-control logic (figs. 1, 2, and 
10) . These equations are nonlinear and are solved numerically in the case of 
the AWJSRA; this is expected to be the case for powered-lift aircraft gener- 
ally so that efficient general approaches to the solution of such equations 
are of interest in the control system design. 

An exhaustive grid search procedure was used in this work and is described 
in this appendix. This method is conceptually simple; data defining the 
aerodynamic and engine forces are assumed given as tables of values on a 
grid of points covering the domain of interest. Table interpolation is used 
to provide a piecewise linear model on each piece of the grid and these linear 
equations are solved piece by piece until the piece containing the solution 
is found or the grid is exhausted. The existence of a solution is determined 
in a finite number of steps, and its uniqueness is determined from properties 
of the model data. 

Alternative methods based on Newtonian iteration were also investigated. 
The classical Newton method is simple and well known but subject to an array 
of difficulties in practice affecting both reliability and computational 
requirements. These difficulties motivate the modifications found in a 
variety of quasi-Newton methods (ref. 13), such as the Newton-Powell algorithm 
(ref. 14). An advantage of iterative searches is their general applicability; 
the same iterative procedure applies to the solution for any pair of variables 
in the AWJSRA trim equation and to solving the trim equation for any powered- 
lift aircraft, including cases where more than two equations in two unknowns 
must be solved, such as the tilt rotor aircraft which requires the solution 
of eight equations in eight unknowns. The Powell algorithm requires little 
more than a user-supplied function evaluation subroutine. In addition, the 
adaptive search of an iterative scheme may converge in fewer steps on a dis- 
tant solution than the rigid grid search pattern of the exhaustive grid search 
method, especially if the aircraft model is not strongly nonlinear. 

For flight control use, important properties of the trim solution algo- 
rithms are computation time and reliability. The Powell algorithm was tested 
with the AWJSRA model after modification to enforce constraints on the search 
domain. It proved efficient with convergence to a solution in a few steps 
for all cases tested. However, methods of proving the convergence of such 
algorithms for all operational conditions, given the aircraft model, do not 
appear to be available. Consequently, the exhaustive grid search procedure 
was adopted for the present control system design. Computation time require- 
ments are proportional to the product of the time required for each step in a 
search and the number of steps to converge on a solution. An upperbound on 
computation time is given by the case of searching for a solution when none 
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exists. This case will occur in practice in the absence of existence condi- 
tions that can be tested before entering the algorithm. When solutions do 
exist, the convergence rate or average number of steps to converge is a prop- 
erty of interest. 

For reliability, convergence of the algorithm to a solution, if it exists, 
and correct selection of a single solution if several exist, are necessary 
properties. In addition, such software properties as insensitivity to round- 
off error, and risk of programming error, which increases with the complexity 
of the algorithm’s structure, are relevant. 

For convenience, the trim equations for the AWJSRA are repeated here 



cos (a + v) 
sin(a + v) 


) 



(Bl) 


where Th(Njj//t, 6, V^) , mE(Nji//r, 6), Ci)(a, 6f, Cj) , Cl(oi, 6f, Cj) are 
scalar functions of the control variables or of auxiliary functions of the 
controls and are given as tabulated data in appendix A. A following 
section of this appendix (Exhaustive Grid Search Solution Method) describes 
exhaustive grid search algorithms for solving equation (Bl) for values of 
the control pairs, {a, v} and {a, 6f}, which yield the commanded applied 
specific forces, (A^^, Aj^^) , given the values of all remaining parameters 
in equation (Bl) . Searches over one and two-dimensional grids are used in 
the solutions. These two algorithms are sufficient for the work described 
in this paper but, in general, solutions of equation (Bl) for any pair of 
variables among {x, p, u} are of interest in the analysis of various 
performance problems. 


Table interpolation is required extensively for function evaluation in 
these algorithms. Interpolation formulas which are continuous over the table 
domain are given next. 


Table Interpolation 

The evaluation of functions, such as Cl, Cq, etc., by interpolation of 
stored tabulated values is used repeatedly in the trim solution algorithms. 
Interpolation formulas for functions of three variables are given below; 
those for functions of one or two variables follow as special cases. 

Assume that a table of values of the function, f(x, y, z) , has been 
stored for all points in a grid covering the domain of interest; that is, 
the set of values 

|f(xi,yj,Zk) ; i = 1, 2, . . . , N^, j = 1, . . Ny, k = 1, . . } (B2) 

This table contains * Ny • Ng values and the grid subdivides the domain 
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^={(x,y,z): xe[x^,XN^], y^[yi»yNy]> 




(B3) 


into (Nx - 1) • (Ny - 1) • (N2 - 1) box-shaped pieces. 

The original function is approximated as piecewise linear by interpolat- 
ing this table. For a given point, (x, y, z) , in the table is interpel- 

lated in two steps as follows. First, locate the grid piece containing the 
point (x, y, z) ; that is, find 


(i,j,k) such that x E y E [yj,yj+i), 2 E (B4) 


This is easily done by several methods; the one used here is defined later 
in equation (B14) and is efficient in real time simulation and control. 
Having found (i, j , k) it is convenient to define the following vector which 
locates the point (x, y, z) relative to y j , z^) and normalizes the 

coordinates 


P = 








(B5) 


The coordinates of p are each numbers in [0, 1). For brevity, the notation 
fiQo> . . ^ 111 ^ adopted to refer to the tabulated values of f 

at the eight corners of the grid piece: (x^, y j , z^) , (x^_j_^ , y j , Z]^) , . . 

Second, the linear approximation of f is given from the usual truncated 
Taylor series expansion of f about (xj^, y j , z]^) : 

" ^000 ^ ^x’'^y’^z ^ 


The approximation, f, is required to be continuous and to match the tabulated 
values- However, the original function, f, is in general nonlinear, and its 
values at the eight corners of the grid piece will not fall in a single 
hyperplane, such^as equation (B6a) . Consequently, to obtain the desired 
properties for f it is necessary to define hyperplanes by calculating the 
gradient vector, Vf, from a different set of four grid points depending on 
the location of p within the grid pieces. This is illustrated in sketch (a) 
for the simpler case of a function of two variables. The rectangular domain 
is divided into two triangular pieces and the gradient calculated from a 
different set of three grid points in each triangle as noted in the sketch. 

The resulting function is linear in each triangle and continuous at their 
common edge. In three dimensions the linear approximation becomes more 
complex and requires subdivision of the box-shaped domain into six tetrahedral 
pieces; approximate formulas for the gradient vector are: 
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Sketch (a) Linear interpolation- 


interpolation formulas other than equations (B6) are also available and 
have different properties; for example, subdivision of the grid piece as in 
equation (B6b) is avoided by adding terms in cross-products of {p^, Py, p^} 
from the Taylor series, which yields 
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+ PxPy^^llO ^010 ^100 ^000^ 
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101 001 100 000 


+ PyPz^^Oll ^001 ^010 ^000^ 


'*' PxPyPz^^lll ^011 ^101 ^001 ^110 '*’ ^010 '*' ^1 


00 


’000 


) 

(B7) 


This formulation extends to more than three variables by retaining all terms 
to first order in {p^^, py, p^, . . .} from the general Taylor series expansion 

(cf . ref . 15) . 

The linear approximation in equations (B6) works well for functions of 
two variables, including the inversion of such functions, but equation (B6b) 
is difficult to generalize to four or more variables. The quasi-linear 
approximation of equation (B7) has worked easily for functions of two or three 
variables, including their inversion. For some applications, it is useful 
for f to have continuous first derivatives over the table domain but 

that property was not required in the present work; equation (B7) has discon- 
tinuous derivatives on the surfaces of the grid pieces while equations (B6) 
has discontinuous derivatives along all interior boundary surfaces implied by 
equation (B6b) as well as on the grid piece surfaces. 


Exhaustive Grid Search Solution Method 


Solution for nozzle and angle of attack— The problem is to solve the trim 
equations (Bl) for {a, v} given the values of all remaining variables in the 
equations. It is convenient to define the quantities 

''' + ™E^a) \ 

^ ' W I 


E 


N 


“ QSw^l 

w 


(B8) 


E E (E^,E^)^ 

i 


In equations (B8) , Cq and a^te functions of a and are assumed tabulated 

on a grid of values 


i 1, • . ., Nqj } 
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where corresponds to stall or maximum The trim equations can be 
rearranged to isolate the nozzle angl^ in a single equation by substituting 
equations (B8) in the trim equation, which becomes 

- T 
E = Tjj(cos(a + v),sin(a + v)) 


and then it follows that 


E(a)l =if 


V = 



(B9) 

(BIO) 


Since E depends on a only and can be tabulated for the grid points, {aj_}, 
equation (B9) can be solved numerically for a after which v is calculated 
from the closed form expression in equation (BIO) . 


The problem of solving equation (B9) is visualized in sketch (b) , which 
shows the operating point (Au^, which trim values of {a, v} are 

desired, and the line traced out by the aerodynamic forces 



Sketch (b).- Solution of equation (B8). 
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as a is varied over the domain [a^, • This line is closely related 

to the Cl - Cj) polars (appendix A) ; it has a near-vertical slope and A]^ is 
strictly increasing in a below the stall angle. The solution of equa- 
tion (B9) is equivalent to finding a* such that the distance from the line 
to the operating point, |E(a*)|, matches the ratio, Th/W; this is visualized 
in the sketch as the intersection of the line with a circle of radius Tjj/W. 

In general, there may be no solution (the line is entirely outside or inside 
the circle) , one solution (one endpoint of the line falls inside the circle) , 
or two solutions as in the case sketched above. When two solutions occur, 
the lower value of a is selected in order that the resulting nozzle angle 
fall within the hardware limits on the nozzle, [6®, 104®]. The higher value 
of a corresponds to a negative value of and this implies a negative 

nozzle angle and downward directed thrust vector from equation (BIO) in 
almost all cases. 


The numerical solution of equation (B9) is considered in two parts: 
first, a method is defined for testing whether a given interval of the a-grid, 
say [a^, o^i+i), contains the solution; and, second, a procedure for searching 
the intervals of the grid is given. 


The interval [a^, ai+i) can be tested as follows. Define the notation 


E. E E(a^) 

AE^« E E . . - E . 

^ i+1 X 

a - a . 


a 


Then E(a) is given on [a^^, by: 


E(a) = + AE^p^ 


► 


J 


(BID 


(B12) 


This is substituted in equation (B9) and the resulting quadratic in Pq^ 
solved : 



While the quadratic has two solutions for p*, only the smaller solution, 
given in equation (B13) , need be tested because the larger one corresponds to 
nozzle settings outside the nozzle hardware limits. If the condition 

pS e [0,1) 

is satisfied by equation (B13) then the solution occurs in the tested inter- 
val, [ot^, is 
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a* = 


a - + 

X 


Pa (“i+i - “1> 


Otherwise, the solution either does not exist or occurs in some other interval 
of the a-grid . 

Each interval of the a-grid is tested as described above until either a 
solution is found or the grid is exhausted with no solution in ” 1 steps. 
While the grid can be searched in any order a minimum computation time search 
ordering is desired for real time control or simulation applications. If a 
solution exists, computation time is proportional to the number of intervals 
tested before finding the solution and can be minimized, on the average, by 
testing intervals in the order of the likelihood that they contain the solu- 
tion. For flight control, solutions are required every control cycle and the 
expected value of the solution corresponds to the nominal path and is either 
constant or slowly varying with time. Consequently, the most probable inter- 
val is the same interval in which it was found for the previous control cycle 
and probability decreases with distance from this interval. The corresponding 
search ordering begins in the same interval, I, in which the solution was 
previously found and expands to neighboring intervals out to the ends of the 
grid; that is, in the order {I ± k} generated as k varies in the expression 


! I-k k=l, ...,I-1 

I + k k = l, ...,No^-I-l 


(B14) 


The maximum computation time required with this search procedure occurs 
when all intervals must be searched and is proportional to - 1. This 

worst-case time is minimized by minimizing the number of grid points in the 
a-grid used in the aerodynamic model within the limits imposed by the desired 
model accuracy. 


Sotut'lon for throttle and angle of attach- The problem is to solve the 
trim equations for {a, 6^} given the values of all remaining variables in 
these equations. The applied acceleration due to engine and aerodynamic 
forces is expressed as a function of (a, by 


a(a,N^) 



cos (a + v) 
sin (a -1 v) 


) 



(B15) 


Here, denotes the corrected engine power, It is convenient to 

solve for Nq in place of throttle since the engine output is tabulated in 
appendix A for this parameter on the grid 

2, . . ., N|^ 

The throttle setting corresponding to any value of is then given from 

figure 30. The problem is now to determine (a*, N*) such that the commanded 
acceleration (A^^, A^^) or, for brevity, a^., is obtained; that is, solve 


62 



a(a,N^,) = 


(B16) 


The applied acceleration, a(a, Nq) , can be tabulated for values of 
(a, N^) in the grids used for the aircraft and engine models. This yields 
the table 



i = 1, 


Nj,, j = 1, . . I (B17) 


On a grid-piece of the domain, say [a^, a ] (g) [Nc j , Ncj_|_;^ ] j a is given 
from the linear interpolation formula (B5) as 


where 


a(a,Nc) = aij + ’ ^ct^PNe ^ 


Pn. 


a - a . 


a - a. 


(B18) 


PNc 




j 


Using equations (B6) , the Jacobian is given by the 2 

L C^i+l,j+l ^i, j+1 j+1 “ 

The problem of solving equation (B16) is visualized i 
tions (B17) and (B18) map values of (a, N^) into values of (A^, A;^) and this 
map is shown in the sketch for the grid lines of the tabulated model. The 
mapping differs depending on configuration; in the conventional configuration 
and for CTOL aircraft generally the effects of the engine power and angle of 
attack are nearly orthogonal and Aj^(a, N^) is nearly independent of engine 
power. Consequently, there is little coupling between the two equations and 
the normal acceleration equation can be solved approximately for a*. However, 
in the powered-lift configuration, Aj^(a, depends strongly on engine power, 

both from thrust vectoring and lift augmentation, and a* can no longer be 
solved separately from nJ, Nevertheless, the solution of equation (B16) , if 
it exists, is unique at all flight conditions and configurations in the design 
domain based on the model properties: 

1. Aj^(a, N^) is monotonic and strictly increasing in a and 

2. A^(a, N(>) is monotonic and strictly increasing in along any 

line for which A|^(a, N^) is constant 


X 2 matrices 


Pn, < Pa 
Pa ^ PNc 


(B19) 


.n sketch (c) . Equa- 
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Sketch (c).“ Acceleration map for (a, 5^). 


The numerical solution of equation (B16) is considered in two parts; 
first, a method of testing whether a given grid piece contains the solution 
is given and then a procedure for searching the pieces of the two-dimensional 
grid is given. 

A grid piece, say [aj_, ot ] (S> [Ncj > Ncj + i ] , can be tested to determine 

if it contains the solution by solving equation (B18) 



(B20) 


This solution is calculated once for each Jacobian given in equation (B19) . 
The result in equation (B20) gives the solution of equation (B16) if it 
satisfies the necessary and sufficient conditions 


0 < 



< 



1 


or 


0 < < 1 

depending on which of the two gradients in equation (B19) was used. If these 
are satisfied, then the solution is 


“* = “i + Pa^^i+l ■ “i) 

Nc = Ncj + - Ncj) 
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Each piece of the grid is tested until either a solution is found or the 
grid is exhausted with no solution in (Nq^ - 1) • - 1) steps. A search 

ordering which is efficient in real-time use is defined next by extension of 
equation (B14) to two dimensions. This ordering tests grid pieces approxi- 
mately in the order of their probability of containing the solution; that is, 
it begins at the grid piece (I, J) in which the solution was found for the 
preceding control cycle, and expands out to the ends of the two-dimensional 
grid, testing groups of grid pieces which frame the starting piece (I, J) at 
each step. The grid pieces in each frame are assumed to be equally likely to 
contain the solution. The indices for the grid pieces in the kth such frame 
(see sketch (d)) are enumerated in equation (B12) for each of the four sides 
of the frame; grid pieces and sides which are outside the model domain are 
deleted from the enumeration by the extremes placed on i, j, and k. 


MODEL DOMAIN AND 


A 

GRIDLINES 







^ k = 0 (GRID PIECE l,J) 
^ k = 1 
k = 2 
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Sketch (d).- Two-dimensional grid search pattern. 




(k) 


(I - k,j) 

(I + k,j) 

' (i,J + k - 1) 
(i,J - k) 


k < I 


k < 


N. 


^ < Na 
k < J 


I 

J 


V 


3 


i 


3 min ’ J min 
■'min 


+ 1 , . 


• »Jmax 


jmin = ~ ; 

jmax = min{NNc - 1,J + k} 

“ ^min’^min 


> -^max 


imin = maxU,I - k + 1} 
imax = mln{N(^ - 1,1 + k - 


1 } 

(B21) 


The search proceeds, taking k = 0, 1, 2, . . ., until a solution is found 
or the grid is exhausted after - 1) • (N]^^ - 1) grid pieces have been 

tested. This search ordering is applicable to the solution of any function 
of two variables. 


I 
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APPENDIX C 


CONTROL MARGIN AND CONFIGURATION SCHEDULE COMPUTATIONS 


Computation of Control Margin 

Control margin is a measure of the acceleration capability available to 
regulate path errors by varying the feedback controls around their trim values 
associated with the nominal flight condition. Control margin (^epends strongly 
on configuration and appears in the configuration schedule design both as a 
constraint on acceptable configurations and as a parameter to lj>e optimized. 
Several different control margins 

CM5j(x,p), CM5^,6t(x»P) . CM6j,ct(x,P) 

were defined in the text associated, respectively, with the usej of three 
controls for regulation and with various pairs of controls. The subscripts 
in the notation indicate the controls held fixed in each case. Control margin 
was defined as follows. The acceleration capability of the regulator, given 
5f, Vg, p, is 

^ 6 f(VE>P) = {(A^(Ve,P,u) ,Ajj(Ve,P,u)) : “ ^ ^ ^ 

where accelerations are generated using the aircraft model (eq.(9)). Next, 
define the region enclosed by ellipses with 5-to-l axes ratios and centered 
at ac 


^(r,Auc,An^) h {(Au.An): 25(Au - Au^)2 + (A^ - A^^)2 < r2 | (C2) 

The control margin is then the largest such ellipse that can be inscribed 
in the envelope of ^6 £(Ve,p), as shown in figure 34. 

CM^jCx.p) E maxjR; (r.A^^.An^) C y^^CVg.p)! (C3) 

Control margin is computed in three steps and these are outlined below. 

1. Generate the envelope of p) (or of 

2. Determine if (A^,, A^^) is inside 5^6 (V£, p) (or inside 

etc.) ^ 

3. If it is, calculate CM 5 ^(x, p) (or etc.) 

The acceleration envelope in step 1 is generated approximately as a 
polygon. The end points of the sides of the polygon are computed as a set of 
points {(A^, A]^)^} on the envelope corresponding to a sequence of values of 

the controls {u^} that form a grid covering the regulator regime. For 
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convenience in later calculations, these points are generated in a clockwise 
manner around the envelope and the controls {uj^} must be taken in a definite 
sequence in order to do so. The required computations are then 


1. Construct grids of values of the controls, a, 6|-, v, covering their 

permitted ranges for regulation, Omax] . t^t^in’ ’ timin’ '^maxl * 

For example , 


a. = a . (a - a . ^ ^ 

1 mxn \ max mxn / n - 1 


6t. = • 

^ X ^mxn 


+ (6 


"•max 


- 6, 


/ n - 1 


"-mxn 


V . = V . + f V 

X mxn \ 

2. Genera t 
following order: 


- V 

max min 


^ ^ 

in / n - 1 


i— 1, . . .,n 


and store 

values of (A^(V£, 

P, 

> An(Ve, 

^i = 

(6f ,a^,6^^,Vj_) 

I— 1 
11 

•H 

. . . , n 

^n+i 


i = 1, 

. . . ,n 

^2n+i 

( 6 £ , , 6 , Vj^) 

i = 1, 

. . . ,n 

^4n+l“i 

(6f , , 6 j v^) 

i = n. 

. . .,1 

^5n+i“i 

(6£,otj^,6 5 

i = n. 

. . .,1 

^en+l-i 

(6^,a^ ,6tj^,v^) 

i = n. 

. . .,1 


The envelope, P) » consists of 6n points { (Au^, An^) , 

i = 1, . • 6n} with successive points connected by straight lines. A 

typical is shown in figure 35. The envelopes, 

are generated similarly except that another control in addition to flap is 
held fixed, and this reduces the number of points to 4n. 


A simple geometric criterion can be used to determine whether the given 
operating point (Au^.> Aj^^) is inside the envelope (step 2 of the computations). 
Let {pj_, i = 1, . . . , 6n} refer to the points generated in step 1 to define 
the figure, and let d± be the angle subtended at the operating point, p^ , 
by the directed line segment joining p^_, (see fig. 36) 


9^ = sin 


-1 


I, ® , 

X x+1 




1+1 


- < ft- < ^ 
2 < 01 < 2 
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I 



where 


1 


Pi - 


■"i+l 



and then 

Pj, e jj/ iff E 0^ = 2ir 
^ ja/ IFF Z 0^ = 0 

Assuming is in sd ^ the control margin can now be calculated. The 

problem of calculating the largest ellipse of the form given in equation (C2) , 
with center at (Au^, ^Nq) which can be inscribed in the envelope, ja/, gener- 
ated by the points { (Au^_, ANi) } is equivalent to finding the radius of the 
largest circle centered at (5Auc» which can be inscribed in the envelope 

, generated by the points {(5Aui, A^j^)}. For simplicity, the notation 
p’ will refer to points (5Au, A^) . The control margin is calculated as the 
smallest distance from the operating point p^ to any one of the sides of 
the polygon as described next. 


Let be the line segment forming the ith side of the envelope. The 

minimum distance, Rj^, from the point, p^, to this segment is the perpendicular 
distance to the segment if the intersection of the perpendicular occurs within 
L^. Otherwise, it is the lesser of the distances to the end points of L^. 

To calculate first define unit vectors, u, n along and normal to the 

line segment Lj_ (see fig- 37). 


pl+i 


■Pi T 

— = (cos a, sin a) 


Ip1+i - Pil 


n = (sin a, -cos a) 

Then the line segment, can be given as 

Li = |p' = + a5; a G [o,|pi+i ” P[\] } 

where a measures the distance along u from the point, p^. The perpendicu- 
lar from p’ intersects the extension of Lj^ at 



“1 = (Pc ■ 

PP 

• 

u 

and then the minimum distance 

to is 




rlpc - pil 

If 

a± 

< 

0 

\ = j l(p;-pl) 

• ii| If 

“1 

G 


1 Ipc - p1+iI 

If 

“i 

> 

- Pil 
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Finally, the control margin is the smallest such distance among the sides of 
the polygon 


CM(p^,) = min{Ri, i = 1, . . • , N} 


Admissible Flap Settings 

The flap schedule was calculated empirically by determining the range of 
acceptable flap settings at each flight condition and selecting the maximum 
such setting. The computational steps required are described below; they 
utilize the previously defined trim solution and control margin algorithms as 
well as the constraint list from the text. 

1. Define grids {Ve^}, {Auj} and {6f^^} which cover the design ranges 
of interest for Vg, A^ and 6 ^ given in 

2. Generate and store acceleration envelopes for these grids 

The standard parameter values, p^, are used. 

3. For each value of x in the set given by 

Ajj = 1 

Ve ^ 

determine those 6f in {6^ } for which 

CM5^(x,Ps) > 0.25 g 

and at least one trim solution exists satisfying the remaining constraints. 
This is done by enumerating solutions for a given throttle grid covering its 
permitted range and evaluating the constraint functions {f^(x, p^, u) } listed 
in the text (fig. 15). 

These steps produced the principal empirical data on which the schedule 
was based; that is, the acceptable range at each point in a grid covering a 
subregion of the design flight regime, , given by 

An = 1, p = Ps 

The flap schedule was selected approximately as the maximum acceptable flap 
at each point and a schedule, which is a function of only two variables, 

A^) , is produced. Step 3 is repeated at other A^ in a grid {Aj^^} 
covering the design range of A^ to determine if the acceptable flight 
envelope for the schedule can be usefully maximized by either adjusting the 
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schedule or including as an independent variable. Similarly, steps 2 

and 3 are repeated over a grid {p^^} covering to determine if p should 

be an independent variable in the flap schedule. 


Optimum Nozzle Schedule 

The optimum nozzle schedule was determined empirically by calculating 
the range of acceptable nozzle values at each flight condition and at the 
scheduled flap setting for that flight condition, and then selecting the 
value with the maximum fixed-nozzle regulator control margin CM^ 5 £^^(x, p) . 

An optimized throttle schedule for maximum fixed-throttle control margin, 
CM(5£^5^(x, p) can be calculated analogously. The required computational steps 
are described below; they utilize the previously defined trim solution and 
control margin algorithms (appendix B) , the constraint list from the text 
(fig. 15), and the flap schedule, F*(Ve, A^) . 

1. Define grids and {v^} which cover the design 

ranges for V;^, A^^, A^ and v '^in 

2. Generate and store fixed-nozzle acceleration envelopes for these grids 

Ps)} 

where the standard parameter values Pg are used and .is taken from the 

the flap schedule, A^^ ) . 

3. For each value of x in the grid, calculate the trim solutions 

corresponding to corresponding values of the constraint para- 

meters {fj_(x, pg, u} and control margin, CM(5£,v(x, Ps)« These data give the 
range of nozzle settings for which trim solutions exist, the subset of this 
range for which they are acceptable, and the acceptable nozzle angle having 
maximum control margin, . 

These steps produced the empirical data from which the schedule, v* (x) , given 
in the text was defined. These data are all taken at the standard parameter 
values Pg, but steps 2 and 3 can be repeated over a grid covering 

to determine if the acceptable flight envelope of the schedule can be usefully 
maximized by including p as an independent variable in the flap schedule. 

The effects of off-nominal parameter values on the flight envelope 
boundaries are significant since the STOL approach attempts to exploit the 
margins of the acceptable flight envelope, particularly 

However, a study of these effects and the possibilities for maximizing the 
acceptable envelope for the designed schedule by including the parameters as 
independent variables was outside the scope of the present effort. 
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APPENDIX D 


ATTITUDE AND ATTITUDE RATE COMMANDS AND STABILITY AXIS 
SPECIFIC FORCE COMMANDS 


The attitude and attitude rate cominands associated with the execution of 
the acceleration commands and the nominal path are calculated in the trim map 
as indicated in figure 25. The Euler angle attitude commands, ($, 0, Y) and 
body axis rate commands (p, q, r) can be calculated from the values of the 
angles {a, 3, their rates for the nominal path and acceleration 

commands, as derived in' this appendix. In addition, the Trimmap requires the 
transformation of the input acceleration commands to stability axis components 
of the applied specific force. This transformation was given in the text 
(eqs. (5)~(7)) assuming zero steady-state sideslip angle. A more general 
formulation for the case of nonzero sideslip is derived in this appendix. 


Attitude Commands 


The notation T/i^g will denote the transformation of vectors referred 
to axes A into vectors referred to axes B. The axes of interest here are iner- 
tial or runway (r) , path (p) , stability (s) , and body (b) . 

The attitude of the aircraft body axes with respect to inertial axes is 
given by the standard Euler angles, 0, '1' , measured by gyros on the aircraft. 
The corresponding transformation is (ref. 10) 


T 


br 


cos 0 cos 

sin $ sin 0 cos f - cos $ sin f 
cos ^ sin Q cos 'F + sin sin 


cos 0 sin '1' 

sin ^ sin 0 sin 4^ + cos ^ cos f 
cos $ sin 0 sin ^ - sin ^ cos 0 


-sin 0 
sin cos 0 
cos ^ cos 0 


(Dl) 


This can also be expressed in the form of a sequence of rotations about a 
single axis (ref. 10) 


= L^(<E.)L2(0)L3(^) 


(D2) 


where the subscript indicates the axis about which the rotation occurs 
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L„(x) 



'1 

0 


0 

= 

0 

cos 

X 

sin X 


-0 - 

sin 

X 

cos X _ 












cos 

X 

0 

-sin X 


0 


1 

0 


_sin 

X 

0 

cos x_ 


cos 

X 

sin 

X 0 “ 


-sin 

X 

cos 

X 0 


0 


0 

1 . 


Expressions for the Euler angles can be given in terms of the angles {a, 3> 
y} by equating equal transformations between inertial and body axes 

T = T T T (D3) 

br bs sp pr 

The transformation to path axes is defined from the air velocity vector 
using equation (6) of the text 

The transformation from stability axes to body axes is simply 

V - 

The remaining transformation from stability to path axes is defined with the 
aid of sketch (e) . The path axes {u, m, 5} are along the air velocity vector 
(u) and normal to it with m and n _in the horizontal and vertical planes, 
respectively. The stability axes {Ig, kg} are^ respectively, along the 

projection of in the body plane of S 3 nnmetry (ib> Jb) > along the body wing 
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axis and in the plane of S3nranetry normal to ig. The stability axes are 

obtained from path axes by first rotating about the air velocity vector direc- 
tion II, through the angle then about the kg axis through the side- 

slip angle 3, in a left-handed direction. 

% = (D6) 

The transformation identity in (D 3 ) can now be given as 

L^($)L 2 ( 0 )L 3 ( 4 ') = L 2 (a)L 3 ’^( 6 )L^( 0 ^)L 2 (Y)L 3 ('Fv) (D 7 ) 


The attitude commands can therefore be calculated from the remaining five 
angles. In the Trimmap, is obtained from the trim solution algorithm, 3 ^ 

is usually taken as zero but can be nonzero in the case of decrabbed flight 
so is retained here for generality, ^vc is one of the cylindrical coordinates 
of the total acc elera tion command, and (y, are given by the commanded air 

velocity vector, VAC. 


To derive formulas for the Euler angles, represent the elements of the 
matrix on the right-hand side of equation (D 7 ) as [^ij ] ^nd then, using 
equation (Dl) obtain 


0 = sin ^-^33) 

7 T ^ 7 T 

2 < 0 < 2 

I* = tan" ^ ^ ^ ) 

/A 

-e 

A 



0 = tan ^ ^ ^ \ 





or, for computational convenience 


(D8) 


$ 


sin 


^23 
cos 0 


The required terms for the calculation of the Euler angles are 

a^2 “ “ ^os a [cos 3 sin y -f sin 3 sin (J)^ cos y] - sin a cos cj)^ cos y 

^23 “ “Siri 3 sin y + cos 3 sin (p^ cos y 

^33 “ a(cos 3 sin y + sin 3 sin (p^ cos y) + cos a cos cp^ cos y 

a^^ = sin 4 :'^[cos a sin 3 cos cp^ - sin a sin <p^] ^ 

+ cos ¥^[cos a (cos 3 cos y - sin 3 sin (p^ sin y) - sin a cos (p^ sin y] 

a^^ = sin ^^[cos a (cos 3 cos y - sin 3 sin <p^ sin y) - sin a cos (p^ sin y] 

+ cos ^^[-cos a sin 3 cos (p + sin a sin (p ] 

V V ^ 

(D 9 ) 
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The heading angle expression can be simplified by denoting the coefficients 
of sin cos ^y in a ^2 as A, B and then we have 


12 

= A 

sin 


+ B 

cos 



11 

= B 

cos 

% 

- A 

sin 


► (DIO) 




^y + tan 


B 

A 




For passenger operations in STOL aircraft, the angles a, Y> © can be 
taken as small angles and second and higher-order terms (H.O.T.’s) in these 
angles can be neglected to obtain the approximation 


0 = Y + a cos + 3 sin + H.O.T.’s 
$ = + H.O.T. 's 

m = a sin - 3 cos + H.O.T. ’s ^ 


(Dll) 


Attitude Rate Commands 


The body axis attitude commands are calculated from the five independent 
angles, {a, 3> Y? ^-\j} > and therefore the body rate commands can be con^ 

structed from the angular rates of these angles and oriented along their 
appropriate axes about which the rotation occurs; that is, from the vector 
identity 


0) = aj^ + 3k^ + + y^t^ 


V r 


(D12) 


All vectors are unit vectors from one or another of the axis frames previously 
defined, and all can be transformed to body axes using equations (D4)-(D6); 
that is, the body axis rate commands are 


^0 


^0 




oV+ (iV 

.0/ Vo/ 


+ ^ T 
V brl 


(D13) 


or 



-sin a 


( cos a cos 3 
sin a cos 


sin 3 


sin a sin (|)^ - cos a sin 3 cos 
cos 3 -cos cj)^ 

-cos a sin - sin a sin 3 cos 



'-cos a[sin y cos 3 + sin 3 sin cos y] “ sin a cos cos y 
-sin Y sin 3 + cos 3 sin (j)^ cos y 
-sin a [sin y cos 3 + sin 3 sin cos y] + cos a cos cf) cos y, 
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The small angle approximation is obtained by neglecting terms which are second 
or higher in y, a, 3, a, 3, Y. ^'v 



The angular rates for y, can be generated from the nominal accelera- 

tion commands, using equations (5) and (7) of the text 







VC 


$ 

VC 


ac * m 
V cos y 

003 ^$ , 

= — (ac • m) + H.O.T. ’s 

g COS y - ac • n dt 


/ 


(D15) 


Of these angular rates has a steady value during turns (of the order of 

3 to 5 deg/sec), varies transiently at turn entry and exit and is zero other- 
wise, and appears in both q, r. The flight-path angle rate y occurs 
transiently in relatively short duration flares and pitchovers and affects q 
principally. The roll rate, occurs only transiently at turn entry and 

exit. It is calculated from the lateral jerk, which is available from the 
command generator as a result of its transition maneuver command generation 
calculations. Of the remaining rates in equation (D14) , 3^^^) assumed 

zero almost everywhere, and can be constructed as the derivative of the 

output of the trim solution algorithm. However, for the nominal path, a 
and y transients are associated but the a transient usually has a higher 
frequency and is less significant as a rate command for following the nominal 
maneuver . 


Stability Axis Specific Force Commands 


The inputs to the basic Trimmap are the stability axis components of the 
applied specific force commands. This applied specific force command, refer- 
enced to path axes, is readily computed from the commanded acceleration ac 
and air velocity vector VAC, using equations (3) and (6) from the text 


( 




A 


m 


A 


n 





(D16) 


This is transformed to stability axis components using (D6) 
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(D17) 



E3(-3)E^(V 



cos 3 “sin 3 cos 
sin 3 cos 3 cos 
0 -sin 


-sin 3 sin 
cos 3 sin 
cos 



These three equations contain five unknowns: 3j ^z* Normally, 

it is assumed that 3, Ay are zero after which A^, A^, can be solved 

with the same result previously given by equation (7). More generally, the 
steady-state value of 3c need not be zero, as in decrabbed flight in the 
presence of lateral winds, and can be assumed given independently of the 
trajectory command. The value of Ay in this case is not arbitrary but 
follows from 3^. and the directional equilibrium of the aircraft; in steady 
state, 3 generates both lateral acceleration (proportional to Cyg3) and 
yawing moment (C^g3) which is balanced by additional lateral force from the 
rudder. The resulting lateral acceleration in steady state is then 


A 


y 




(D18) 


where b, are, respectively, the wingspan and the moment arm of the 

rudder center of pressure about the aircraft center of gravity. Consequently, 
equation (D17) can now be solved for assuming that A^, A^, A^ 

and 3 q are given and that Ay is computed from 3^. The solution is 


where 


A - sin 3 A_r 
^ _u c T 

^ cos 3c 

^z = - I 



E 


A 

y 


- sin 3^A^ 
cos 3^. 


sin ^ 


A* 


+ An 


2 


(D19) 


Equations (D16) , (D18) , and (Dl9) comprise the map of the applied specific 
force commands to stability axes used in the Trimmap of figure 25. 
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TABLE 1.- TABULATED AWJSRA ENGINE MODEL - ONE ENGINE, 
CORRECTED HOT THRUST, COLD THRUST, INLET MASS FLOW 
AND FUEL FLOW AND THROTTLE - POWER RELATION 



0 

53.6 

84.0 

89.5 

92.5 

95.0 

98.0 

103.515 

Tabulated 

parameter 

Knots 










0 

0 

380 

4552 

10,688 

16,521 

21,506 

26,098 

31,198 


60 

0 

509 

4858 

10,928 

16,740 

21,692 

26,256 

31,325 

% 

^ “V N 

250 

0 

2473 

9350 

14,368 

19,911 

24,309 

28,518 

33,125 

0 

0 

0 

363 

3444 

7385 

11,096 

13,246 

14,684 

15,712 


60 

0 

378 

3486 

7414 

11,074 

13,142 

14,632 

15,809 

_ ^ N 

250 

0 

596 

4068 

8039 

11,870 

13,836 

1 

15,315 

16,380 

0 


■ ' 1 

0 

19.8 

36.7 

58.4 

72.4 

- 

80.4 

■ ■ i 

86.6 

'1 

94.7 

nig kg/ sec 


0 

551.6 

864.7 

1401.0 

2206.9 

2452.7 

2965.0 

3782.9 

-E 

^ kg/hr 

6/r” 


Note: Throttle-power relation: 

= max(0.0618 N^, 1.3866 N^ - 107.17) 
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TABLE 2 


AWJSRA DRAG AND LIFT COEFFICIENTS 


( a ) C ^ 


a 

- 10.5 

- 6.5 

- 2.5 

1.5 

5.5 

9.5 

13.5 

17.5 

19.5 

27.5 






Flap 

= 5.6 





0 

0.092 

0.068 

0.073 

0.095 

0.136 

0.188 

0.261 

0.375 

0.443 

0.811 

.2 

-.111 

-.135 

-.124 

-.105 

-.062 

-.005 

.070 

.200 

.269 

.648 

.4 

-.313 

-.335 

-.324 

-.299 

-.256 

-.198 

-.117 

.018 

.092 

.484 

.6 

-.510 

-.532 

-.520 

-.495 

-.446 

-.378 

-.297 

-.161 

-.087 

.321 

.8 

-.708 

-.730 

-.718 

-.688 

-.638 

-.570 

-.478 

-.337 

-.257 

.146 

1.2 

- 1.101 

- 1.127 

- 1.116 

- 1.078 

- 1.023 

-.949 

-.846 

-.690 

-.610 

-.184 

2.0 

- 1.883 

- 1.905 

- 1.901 

- 1.850 

- 1.780 

- 1.674 

- 1.562 

- i .:^94 

- 1.304 

-.847 






Flap 

= 30.0 





0 

0.040 

0.110 

0.160 

0.220 

0.270 

0.330 

0.370 

0.410 

0.430 

0.500 

.2 

-.240 

-.150 

-.070 

.020 

.110 

.190 

.270 

.340 

.380 

.510 

.4 

-.480 

-.370 

-.270 

-.160 

-.060 

.040 

.140 

.240 

.290 

.480 

.6 

-.700 

-.580 

-.470 

-.350 

-.230 
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Figure 2.- AWJSRA Force Trimmap 
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(a) Transformation angles. (b) Path and runway axes. 

Figure 3.- Axis systems and transformation angles. 
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Figure 4.- Applied specific force commands - steady flight conditions 




Figure 5.- Applied specific force commands - transient commands for flight- 

path angle changes. 
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Figure 6.- STOL test approach path: hodograph of applied specific force 

command . 
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Figure 7.- Augmentor wing jet STOL research aircraft (AWJSRA) . 
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(b) Glide slope. 


Figure 9.- Concluded 
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(b) Basic Trinraiap - AWJSRA. 

Figure 10.- Basic Trimmap structure - AWJSRA. 
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Figure 11.- Regulator control regime - fixed nozzle mode; 
Vg = 65 knots, 6^ = 65°, v = 84.2°, p = p^. 
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Figure 14.- Flow diagram - Trliranap regulator control logic 
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(a) Constraint list . 



(b) Flap extremes. 

Figure 15.- Operational constraints on control settings - AWJSRA 
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(e) Upper limits on angle of attack for lift margin 
constraint and regulator usage. 

Figure 15.- Concluded. 
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(g) = -0.18 



Vg, knots 


(d) = 0.15 (h) F*(A^, Vg) 


Figure 16.- Admissible flap range - standard case (p 
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Figure 17.- Nominal flap schedule; contour plot 
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Figure 18.- Flap command rates. 




Figure 19.- Control regimes for fixed-nozzle mode; Vg = 65 knots, 

(Sf = 65“, p = Ps- 
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generating functions at = 1. 


Figure 22,- Nozzle schedule for maximum control margin. 
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(b) Nozzle schedule - contour plot (A 


Figure 22.- Continued. 
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(b) Angle of attack, deg. 


Figure 23.- Continued 
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(c) Trim elevator setting, deg. 


Figure 23.- Continued. 
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Figure 24.- Trim solutions for the configuration schedule - 
Vg = 80 knots (p = Pg) . 
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(b) Nozzle schedule contour plot (Vg = 80 knots) . 
Figure 24.- Continued. 






TRIM SOLUTION 
EXISTENCE BOUNDARY 





(e) Trim elevator setting contour plot (V^ = 80 knots) . 
Figure 24.- Concluded. 
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Figure 25.- AWJSRA Force Trimmap. 
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Figure 26.- Four-dimensional simulation test approach path. 
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(a) Acceleration and speed commands 
Figure 27.- Simulation test path. 
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(c) Trajectory errors 
Figure 27,- Concluded 
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V = AIRCRAFT VELOCITY VECTOR 
= STABILITY AXES 

a.p = ANGLES OF V,T VECTORS FROM BODY LONGITUDINAL AXIS 
L, D = WING-BODY AERODYNAMIC FORCE COMPONENTS 
Lj = TAIL LIFT 

'''h = ENGINE RAM DRAG AND THRUST FORCES 

Figure 29.- Engine and aerodynamic forces - AWJSRA. 
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30.- AWJSRA engine model - one engine, 
mass flow, fuel flow and throttle 


hot thrust , cold 
-power relation. 


thrust. 
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(a) Effect of altitude 
and airspeed on 
maximum continuous 
hot thrust. 

Figure 31.- Engine output 


(b) Effect of tempera- 
ture on maximum 
continuous hot 
thrust. 

at maximum continuous power 
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(c) Effect of air- 
speed on Cj at 
maximum continuous 
power . 

(N^ = 95.8): two 
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a CONTOUR VALUES: -10.5°, -6.5°, -2.5°, 1.5°, 5.5°, 9.5°, 13.5°, 
17.5°, 19.5°, 27.5° 

Cj CONTOUR VALUES: 0, .2, .4, .6, .8, 1.2, 2 

Figure 32.- Lift-drag polars — AWJSRA. 
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Figure 35.- Acceleration envelope, (Vg, p) ; <5f = 65°, Vg = 65 knots 



Figure 36.- Test for existence of control margin. 



Figure 37.- Minimum distance from p’ to a line segment 




1. Report No. 

NASA TP-1222 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

CONFIGURATION MANAGEMENT AND AUTOMATIC CONTROL OF AN AUGMENTOR 
WING AIRCRAFT WITH VECTORED' THRUST 

5. Report Date 

March 1979 

6. Performing Organization Code 

7. Author(s) 

Luigi S, Cicolani, B. Sridhar,* and George Meyer 

8. Performing Organization Report No. 
A-7099 

10. Work Unit No. 
505-07-11 

9. Performing Organization Name and Address 

NASA Ames Research Center 
Moffett Field, California 94035 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Paper 

12. Spojworing Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 

*NRC Postdoctoral Research Associate 


16. Abstract 

An advanced structure for automatic flight control logic for powered-lift aircraft operating i 
terminal areas is under investigation at Ames Research Center. This structure is based on acceler^ 
tion control; acceleration commands are constructed as the sum of acceleration on the reference 
trajectory and a corrective feedback acceleration to regulate path tracking errors. The central 
element of the structure, termed a Trimmap, uses a model of the aircraft aerodynamic and engine 
forces to calculate the control settings required to generate the acceleration commands. 


This report describes the design criteria for the Trimmap and derives a Trimmap for Ames exper 
imental augmentor wing jet STOL research aircraft. The principal problems are associated with con- 
trol redundancy (there are two more controls than necessary to generate any given acceleration 
command) and model nonlinearity. Control redundancy is resolved using a stored configuration sched 
ule which selects two of the controls (flap and engine exhaust nozzle) as a function of the referen 
flight condition while the remaining two controls (elevator and throttle) are computed on the basis 
of the total acceleration command and the aircraft force model. The configuration schedule is 
derived for maximum control margins. The algebraic nonlinearity of the force model is treated usin 
piecewise linear descriptions of the model over its bounded domain. 


The automatic control system, including the Trimmap described in this paper, was subjected to 
simulation tests using a rigorous STOL approach trajectory. These tests demonstrate the system 
response to maneuver commands, Trimmap model errors, and steady winds. The proposed Trimmap 
maintains an optimal configuration and coordinates all controls during any admissible maneuvers 
and in response to steady winds; it compensates automatically for model errors. 


17 . 


Key Words (Suggested by Author (s) ) 
Automatic V/STOL flight control 
Automatic configuration schedule 
Trimmap 

Powered-lift aircraft 


18. Distribution Statement 
Unlimited 


STAR Category - 08 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 
141 


22. Price* 
$ 6.00 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA~Langley , 


THIRD-CLASS BULK RATE 


National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



1 1 10, A, 021779 S00903DS 

DEPT OF THE AIE FORCE 
AF HSAPOHS LABOEATOfil 
ATTN: TICH8ICAL LIBEABI (SOL) 
KIBTLAND AFB 87117 




POSTMASTER: 


If Undeliverable (Section • 
Postal Manual) Do Not Re 


